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1 

INTRODUCTION 

This report describes the present status of several closely re- 
lated, ongoing investigations of the use of synthetic aperture radar 
(SAR) data for oceanographic purposes. The common element in these 
investigations is that they make use of information which is con- 
tained within the data recorded by the SAR system but which is not, 
for the most part, preserved in the processed SAR imagery. The ob- 
ject of these investigations is to obtain information on the motion 
of the water surface, including currents, orbital velocities and 
phase velocities of surface gravity waves. Data analysis methods 
and results are presented for various SAR systems, including both 
aircraft and satellite platforms. 

Two general types of analysis are described in this report. 
Section 2 provides the background for both analysis techniques by 
discussing the general principles of SAR data collection and proc- 
essing. Section 3 describes the Doppler analysis methods which are 
being investigated for the extraction of current and wave orbital 
velocity information. Section 4 describes the use of Doppler filter- 
ing techniques for measuring phase velocities. Section 5 summarizes 
the present status of these investigations and makes recommendations 
for further research. 
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BACKGROUND 

A brief description of a synthetic aperture radar (SAR) system, 
including the data collection and processing methods used for recon- 
structing an image of the scene, is presented in this section. A 
more detailed description of the processing methods for extracting 
current and wave velocity information is contained in the following 
sections. The system described here uses optical recording and proc- 
essing methods, but the basic principles involved apply to both op- 
tical and digital systems. 

A SAR system consists of a pulsed microwave transmitter and re- 
ceiver mounted in an airborne or spaceborne platform, a recording 
system, and a ground-based processor. The transmitter sends out a 
series of short coherent pulses in the broadside direction. These 
pulses are scattered from the surface and the reflected pulses are 
received and recorded by the SAR. By time-gating the received 
pulses, using pulse compression techniques in some cases, a fine 
resolution is obtained in the cross-track direction. 

The method by which the scene is resolved in the along-track 
direction is the principal difference between synthetic-aperture and 
real-aperture radars. In a real-aperture radar, each pulse is taken 
as an independent sample of the scene. The along-track resolution 
is therefore determined by the antenna beamwidth, which is inversely 
proportional to the physical length of the antenna. Thus, a large 
antenna is needed to achieve an acceptable resolution in the along- 
tracK direction. In a synthetic aperture radar, a comparable resolu- 
tion can be obtained using a much smaller antenna by using informa- 
tion from successive pulses to reconstruct an image of the scene. 
This process, which may be viewed in analogy with either pulse com- 
pression or holographic techniques, results in an along-track resolu- 
tion which theoretically approaches one half of the physical antenna 
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length [Brown, 19671. As will be shown in the following sections, 
the fact that a given object is illuminated by a large number of 
successive pulses rather than a single pulse also allows information 
on the Velocity of the scatterers to be extracted from the SAR data. 


The process by which the SAR data is recorded and processed may 
be carried out using either digital or optical (analog) techniques. 
Optical processing techniques are attractive because of the high 
density of the storage medium (film) and high processing rates. The 
advantages of digital processing are the greater dynamic range and 
fidelity of the storage medium and the fact that the processing can 
be more tightly controlled to avoid noise and nonlinearities. 


A schematic diagram of a typical optical processor is shown in 
Figure 1. A photographic record of the raw data received by the SAR 
is placed at the signal plane (p.j) and illuminated with a coherent 
light source. As recorded on film, the signal history for a point 
object forms a Fresnel zone plate which focuses the light on a point 
in the image plane (P 2 ). The purpose of the optical elements be- 
tween the signal plane and the image plane is to bring the light from 
all parts of the signal film into focus at a common plane [Kozma, 
et al., 1972; Cindrich, et al . , 1977]. The light intensity distribu- 
tion at the primary frequency plane (Pp is essentially the 
Fourier transform of the signal plane, thus, by placing a filter or 
aperture in this plane, as discussed in Section 4, a portion of each 
Doppler history can, in effect, be selected for processing. 

For the purpose of generating an image of the scene, a photo- 
graphic film is normally placed at the image plane. However, for 
the measurements described in Section 3, the film is replaced by an 
aperture which defines the sub-area of the scene selected for Doppler 
analysis. Light falling on this aperture is passed through another 
spherical lens (f^) to the secondary frequency plane (P^) where 
its intensity distribution is measured and recorded. 
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As will b3 shown in Section 3, the distribution of the light in- 
tensity at ?2 in the along-track direction contains information on 
the radial (i.e., line-of-sight) velocity of the scatterers within 
the image-plane aperture. The shape of this distribution is deter- 
mined primarily by the antenna gain pattern, and the distribution is 
shifted by an amount proportional to the radial velocity of the scat- 
terers. Thus, if the amount of this shift can be measured, the 
velocity of the scatterers can be inferred. 


6 


1 : 


RADAR AND OPTICS DIVISION 


3 

CURRENT AND ORBITAL VELOCITY MEASUREMENTS 

The sensing of surface currents using conventional SAR data is 
potentially feasible using one of three techniques. One possible 
technique is to use the azimuthal displacement in the SAR image of 
an object which has a velocity component in the range direction 
[Raney, 1971]. The use of this phenomenon for measuring ocean cur- 
rents was proposed by Shuchman, et al . [1979] although no data was 
presented. Measurements of the current in the Columbia River 
(Oregon) was reported by Shemdin, et al . [1980] using this technique 
with Seasat data. By measuring the azimuthal shift in the apparent 
boundary of the current as it flowed past the jetties, the radial 
line-of-sight velocity (V^) of the Columbia River was successfully 
estimated. This method may be feasible where sharp current bound- 
aries exist, but is probably not useful in the open ocean where cur- 
rent boundaries are more diffuse. 

A second technique to measure currents from SAR data is summar- 
ized by Hayes and Shuchman [1981] , who showed it is possible to 
calculate the magnitude and direction of the Gulf Stream by observing 
gravity wave refraction across the Stream. This method appears to 
be feasible in deep-water areas where a suitable gravity wave field 
exists, but fails where such a wave field does not exist or where 
wave interactions occur with the bottom. 

The third technique to measure currents reported by Shuchman, et 
al. [1979]* and Gonzalez, et a1. [1981]*, and the subject of 
this report, utilizes a measurement of motion-induced Doppler pertur- 
bations in the SAR signal history. As described in the following 
section, this technique requires that the radial (line-of-sight) 

*Thes e”Tr t fcT^ supported by this NOAA contract are found in the 
appendices to this report. 
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component of the current be large enough to produce a measurable 
shift in the Doppler spectrum of the SAR signals. A similar method 
may be used to infer the orbital velocities of gravity waves by mea- 
suring the broadening of the Doppler spectrum. j 

i 

Data from both aircraft and satellite (i.e., Seasat) systems has [ 

i 

been analyzed using Doppler techniques. These data sets are described 

in Section 3.2, and the results of the analyses are presented in [ 

1 

Sections 3.3 and 3.4. Finally, sources of error in these measure- 
ments are discussed in Section 3.5. 

f- 

3.1 THEORY OF DOPPLER ANALYSIS METHODS 

The theoretical basis for the Doppler measurement techniques can 
be understood by considering the simple case of an isolated point re- 
flector. The geometry of this situation is illustrated in Figure 2. 

The object is assumed to enter the antenna illumination pattern at a 
time t = -T/2 and is broadside to the SAR platform at t = 0. The 
slant range from the antenna to the object is given, to second order 
in t, by 

R(t) = Ro - V r [<'' - <' > 

where R^ is the range at t = 0, V is the platform velocity, 
is the radial velocity of the object, is the along-track veloc- 
ity of the object, and a^ is the radial acceleration of the ob- 
ject. The effects of and a^ are neglected in the following 
discussion, leading to the further simplification 

R(t) * R„ - V^t + (2) 

0 

Note that the variation in R(t) over the imaging time T is small 
compared with the swath width but is large compared to the radar 
wavelength. The SAR transmits a series of pulses at time intervals 
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t = nAt, receives the reflected pulses, and compares the phase of 
the received signals with a reference signal of the same frequency. 
The resulting signal for a point object of unit cross section can be 
expressed as 



Si(t) 



exp [j2kR(t)] 


(3) 


where k = 2tt/x is the radar wavenumber and A(^) is proportional to 
the antenna gain pattern in the along-track direction. In the case 
of optical recording, the real part of this signal is recorded on 
film at the position (x', y') where x' corresponds to the time t and 
y' corresponds to the across-track coordinate. 

The digital or optical processing of the SAR signal in order to 
reconstruct an image of the scene may be described in the ideal case 
as a convolution of this signal with the function 

q(t) = (4) 

P 

where b is normally chosen as b = kV /R^ to cancel out the quad- 
ratic phase factor in S-|(t). In the case of optical processing, 
this operation is carried out by simply illuminating the signal film 
with coherent light and measuring the light intensity at an appro- 
priate distance from the signal film [Kozma, et al., 1972; Cindrich, 
et al., 1977]. The illumination is assumed for the present to be 
uniform and the illuminated area is assumed to be much longer than 
the signal history length. The effects of variations in the illumi- 
nation and of finite aperture widths are considered in Section 3.5. 
The processed signal (corresponding to the distribution of the light 
amplitude in the along-track direction at the image plane of the op- 
tical processor) is then given by 
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00 


S,(t) 




)q(t' - t) dt' 


-jbt^ f /Vt'\ j2(bt-kV )t' 

J it)‘ 


dt' 


2 

c A[2k(Vt - aX)] 


( 5 ) 


where A denotes the Fourier transform of A, c is a complex constant, 
and 



( 6 ) 


Thus, the image of the point object is essentially the Fourier trans- 
form of the antenna pattern displaced by a distance &X which is pro- 
portional to the radial velocity of the object. This displacement 
could be used to infer the radial velocity of the object if the 
proper location of the object were known. However, in most cases, 
this method is not feasible for mapping ocean currents since the 
current boundaries are not sufficiently sharp and are accompanied by 
surface roughness changes which tend to obscure any image displace- 
ment. A more promising approach is to consider the Fourier transform 
of the output signal 

00 

= j Sglt) e^“^ dt . (7) 

—00 

Using the above expression for S2(t), it can be shown that 

= c'|F(» - 2kV^)|^ (8) 

where F(u>) is the Fourier transform of the function 

f(t) . a(«) . (9) 
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In the optical processor, lS 2 (to)| corresponds to the light inten- 
sity in the secondary frequency plane (P^)* with the u) axis lying 
parallel to the longitudinal axis of the film (i.e., the along-track 
direction). Thus, the spectrum observed in the frequency plane is 
shifted by an amount 2kV^ rad/sec, or 2V^/x Hz, for an object 
having a radial velocity V^. If this spectrum shift is an appre- 
ciable fraction of the spectrum width, it can, in principle, be de- 
tected and used to infer the radial velocity V^. , 


The shape of the spectrum is determined primarily by the antenna 
gain pattern in the along-track direction. If the antenna is a uni- 
formly illuminated linear aperture of length D, the one-way (inten- 
sity) gain pattern is given [Brown, 1967] by 



( 10 ) 


where p = \/D is the antenna beamwidth (the 4 dB points are at ^ = 
±S/2). If this function is approximated by a Gaussian 


40 ^ 

" 7 " 

g(0) = e ® 


the spectrum has the form 


|F{(o)j = c e 


- 4^/0 


( 11 ) 


( 12 ) 


i.e., the spectrum also has a Gaussian form with an equivalent width 


0 ' 



Pi 1/2 

+ (kV0)^ ss 2kV0 


(13) 


at the 4 dB points. 

The current detectability criterion using a Gaussian antenna 
pattern can be stated as 
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frequency shift _ 

spectrum width 3v ^ ' 

where a is a threshold whose value is determined by the system per- 
formance and the size of the sample used for measuring the spectrum 
(both of which influence the amount of "noise" in the spectrum mea- 
surement). The data presented in the following sections indicate 
that values of a ~ 0.1 are achievable for areas on the order of 
10^ resolution elements using optical methods. Preliminary expe- 
rience with all-digital processing suggests that at least an order 
of magnitude improvement in the detectability threshold can be ob- 
tained with digital methods. 

The foregoing discussion has assumed that all the scatterers 
within the test area are moving with a constant velocity, as in the 
case of a uniform surface current. When the test area contains sur- 
face elements moving at different velocities, the resulting spectrum 
may be described as a convolution of the uniform-velocity spectrum 
with the velocity distribution function for the surface. This would 
be the case, for example, when the test area contains a gravity wave 
field, because of the orbital motions associated with these waves. 
If the radar cross section were constant over the entire wave sur- 
face, the orbital velocity would cause a symmetrical broadening of 
the spectrum. In fact, however, the radar cross section is modulated 
by variations in both the slope and the surface roughness along the 
wave profile. Describing this effect by the modulation index m 
[Keller and Wright, 1975], the radar cross section at a given point 
on a wave of frequency may be written as 

<j{x, y, t) = 0^[1 + m cos (Kx cos + Ky sin o^ - fit + i^)] (15) 

2 

where is the average cross section, K = fi /g is the water 

wave number, & is the wave propagation direction relative to the 
w 

x-axis, and <j) is the phase angle by which the maximum cross section 
leads the wave crest. The orbital velocity at this point on the wave 
surface is (to first order in the wave amplitude) 
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V„(x, y, t) = a cos e cos (Kx cos a + Ky sin 9 - f2t) (16) 


in the cross-track direction, and 

V,(x, y, t) = a sin (Kx cos 0 + Ky sin 0 , . - nt) (17) 
^ w w 

in the vertical direction. Thus, the radial component of the orbital 
velocity for an incidence angle 0 ^. is 

V^(x, y, t) = sin 0 . + cos 0 . . (18) 

The wave-broadened Doppler spectrum can be obtained by integrating 
over one complete wave cycle, either in time or space. Integrating 
over time, we obtain 


„/ F[“ - dt (19) 

where the position dependence has been deleted from and o, 
since it is cancelled out by the integration process. For a Gaussian 
antenna pattern, the spectrum may be written as 


J 'll cos (nt + 0)] dt 

^ ^ ( 20 ) 
This spectrum is plotted in Figure 3 for three values of 2kJ2a/e' with 
= 0 ^ = 45°, m = 0.5, and ^ = 45°. Note that a significant 
broadening is evident for the case 2kfia/e' = 0.5. The only effect 
of the modulation (m) appears to be a slight asymmetry and shifting 
of the spectrum. 

As in the case of Doppler shifting, the measurable increase in 
the Doppler bandwidth due to wave orbital motions depends on the 
system performance and the sample size. Using = 2kpV, we may de- 
fine the minimum measurable wave amplitude as that which satisfies 


2ir/fi 

2 ^ J 
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FIGURE 3. EXAMPLE OF DOPPLER SPECTRUM 'BROADENING DUE TO WAVE ORBITAL 
VELOCITY, FOR THREE VALUES OF THE SCALED WAVE AMPLITUDE 

2kfia/3- « 

(Frequency axis is scaled by the stationary-surface 
Doppler bandwidth S''.) 
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where a' Is a broadening threshold whose value is on the order of 
0.5. Examples of wave-broadened Doppler spectra are given in Sec- 
tion 3.3 below. 

3.2 DATA SET DESCRIPTION 

Doppler measurements were made on the aircraft and Seasat data 
sets listed in Tables 1 and 2. These data sets are described briefly 
in this section, and the measurements themselvijs are discussed in 
the following section. 

The first four data sets listed in Table 1 were collected with 
the ERIM X-L band SAR system mounted in a C-46 aircraft platform. 
The first three data sets include an area just off the coast of 
Florida, while the fourth covers an area further offshore which in- 
cludes the western boundary of the Gulf Stream. Surface conditions 
for the first three data sets are characterized by a dominant wave- 
length of 80 m with a significant wave height of 1.56 m, and an on- 
shore current of 85 cm/sec. The direction of propagation of the 80 
m waves was also in the westerly or onshore direction. The Gulf 
Stream current velocity distribution on 15 December was such that 
the difference in current velocity for the two areas selected from 
pass 13 was on the order of 0.5 m/sec. 

The remaining data sets listed in Table 1 were collected with 
the same radar system mounted in a Convair 580 aircraft. These data 
sets include two inlets into the Straits of Georgia east of Vancouver 
Island, British Columbia. The 27 July pass was collected during 
flood tide, at which time a current of approximately 2 m/sec was 
flowing into the Straits of Georgia. At the time of the 28 July 
overpass, the current was in the opposite direction at about 1 m/sec. 
The sea state was fairly calm in these areas on both days, with no 
visible gravity waves present in the radar imagery. 

The Seasat data sets considered in this study are .listed in 
Table 2. The Columbia River pass was collected at ebb tide, with a 
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TABLE 2 

SUMMARY OF SEASAT DATA SETS SELECTED FOR DOPPLER ANALYSIS 


Area 

Date 

Orbit 

Look Direction 

Columbia River 

7 July 1978 

150 

66“ 

Dover Straits 

19 August 1978 

762 

63“ 

Gulf Stream 

n August 1978 

651 

63“ 

Straits of Juan de Fuca 

13 August 1978 

681 

63“ 
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current of 2. 5-3.0 m/sec flowing out of the river. The SDR data for 
this pass shows an unusually large yaw angle of approximately 3.5 
degrees. The current in the English Channel (Dover Straits) at the 
time of the 19 August overpass was estimated to be 0.5 to 2.0 m/sec 
at 210’N, based on tide tables for this area. 

The Gulf Stream current on 11 August 1978 is not precisely known, 
but the USCG Weekly Sea Current chart for 30 August 1978 indicates a 
maximum current of about 2 m/sec in an approximately northward direc- 
tion for the area covered by Orbit 651. Striations on the Seasat 
imagery indicate that the Gulf Stream appears to be flowing in a 
somewhat more easterly direction, at an angle of about 45" to the 
range direction. 

For the Straits of Juan de Fuca, tide tables indicate a current 
of approximately 0.5 m/sec into the Straits. 

The incidence angle for all of the Seasat data is approximately 
20", so that the radial component of the current is about one third 
of the horizontal current in the cross-track direction. 

3.3 CURRENT MEASUREMENTS 

Measurements of the Doppler spectrum for selected sub-areas of 
each data set were made on the ERIM optical processor described in 
Section 2. Sub-areas were selected for each data set by placing a 
rectangular aperture in the image plane of the processor, and the 
distribution of the light intensity in the along-track direction of 
the frequency plane was measured and recorded. The aperture size 
was equivalent to a 1 km x 1 km area for the 1975 aircraft data sets. 
For the 1978 aircraft data sets, the aperture size was 250 meters in 
azimuth by 500 meters in range. For the Seasat data, the aperture 
size was typically 3 km x 3 km. 

The frequency scale for these plots was established by inserting 
a grating with a known line spacing at the signal plane and measuring 
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the separation of the resulting spikes in the frequency plane. This 
spatial frequency was then converted into a temporal frequency by 
multiplying by the film speed. The film speeds for the aircraft data 
sets are shown on Table 1. For Seasat, the film speed was 3.24 
cm/sec. 

The Doppler spectra for two sub-areas of the Marineland pass 1 
X-band (3.2 cm) data over water and land are shown in Figure 4. Al- 
though there is a considerable difference in the width of these spec- 
tra, the relative shift between the water and land spectra is neglig- 
ible. The spectrum width, defined as the difference between the 3 
dB points, is the 72.8 Hz for the land spectrum and 160.8 Hz for the 
water spectrum. The L-band (23.8 cm) spectra for water and land are 
shown in Figure 5. These spectra are nearly the same in both band- 
width and mean frequency. 

The X-band spectra for passes 3 and 7, which were flown at ±90° 
with respect to the pass 1 flight path, are shown in Figures 6 and 7. 
The relative shift between the water and land spectra is 45.8 Hz for 
pass 3, and -54.3 Hz for pass 7. Using the equation 

2Vr 

Af = -/ (22) 

these Doppler shifts imply radial velocities of 0.73 m/sec and -0.87 
m/sec, respectively. Dividing by sin 60°, one obtains a horizontal 
current of 0.85 m/sec for pass 3 and 1.00 m/sec for pass 7. The 
direction of this current is indicated by the sign of the Doppler 
shifts for passes 3 and 7, and by the absence of any shift for pass 
1, to be in the westerly direction. 

The L-band spectra for pass 7 (Figure 8) show a negligible 
Doppler shift, despite the existence of an apparent radial velocity 
component as indicated by the X-band data. However, if one calcu- 
lates the expected Doppler shift at L-band due to this velocity com- 
ponent, one obtains a result of 6 Hz, which is apparently too small 
in comparison with the bandwidth to be measurable. 
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FIGURE 4. DOPPLER SPECTRA FOR MARINELAND (PASS 1) X-BAND DATA. 
(Dashed line indicates land spectrum.) 
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FIGURE 5. DOPPLER SPECTRA FOR MARINELAND (PASS 1) L-BAND DATA. 
(Dashed line Indicates land spectrum.) 
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FIGURE 6. DOPPLER SPECTRA FOR MARINELAND (PASS 3) X-BAND DATA. 
(Dashed line indicates land spectrum.) 


23 


Intensity 





Frequency (Hz) 


FIGURE 7. DOPPLER SPECTRA FOR MARINELAND (PASS 7) X-BAND DATA. 
(Dashed line indicates land spectrum.) 
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FIGURE 8. DOPPLER SPECTRA FOR MARINELANO (PASS 7) L-BANO DATA. 
(Dashed line indicates land spectrum.) 
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Doppler spectra for two sub-areas of the Gulf Stream data set 
are shown in Figure 9. These two sub-areas are spaced approximately 
3.1 km apart, and are on opposite sides of an apparent current bound- 
ary in the image. The relative shift between these spectra is 23^4 
Hz, which indicates a difference in radial velocity of 0.37 m/sec or 
a horizontal current of 0.53 m/sec toward the north for the second 
sub-area inside the current boundary. This result is consistent with 
the known direction of the Gulf Stream and with surface observations 
of the current gradient at this location. 

Figure 10 shows two Doppler scans for the Porlier Pass data sets, 
one from a sub-area within the channel itself and one from an adja- 
cent land area. The spectrum width for this data is larger than that 
for the 1975 data sets because of the higher aircraft speed of the 
Convair 580 as compared with the C-46. This increase in the band- 
width, along with the lower signal levels due to the smaller size of 
the sub-areas, causes a noticeable asymmetry in the recorded spectra. 
This effect was modeled as an additive noise component which varies 
linearly with frequency. A straight line was drawn between the high- 
frequency asymptote and the zero-frequency background, and points 
along this line were graphically subtracted from the plotted curves. 
The half-power points were then located on the "corrected" curves 
and used to calculate the bandwidths and mean frequencies of the 
spectra. The Doppler shift of 36.4 Hz between the land and water 
spectra in Figure 10 implies a radial velocity of 0.58 m/sec in the 
flood direction (i.e., into the Straits of Georgia). Assuming this 
current to be at an angle of 45° to the cross-track direction, the 
calculated current speed is 1.43 m/sec. This is smaller than the 
maximum reported current of about 2 m/sec, possibly because it repre- 
sents an average value over an area containing currents lower than 
the maximum. 

The Doppler scans for the three Active Pass data sets on 28 July 
1978 are shown in Figures 11-13. Passes 10 and 11 are at the same 
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. DOPPLER SPECTRA FOR GULF STREAM (PASS 13) X-BAND DATA. 
(Dashed line indicates stationary water spectrum.) 
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FIGURE n. DOPPLER SPECTRA FOR ACTIVE PASS (PASS 10) X-BAND DATA. 
(Dashed line indicates land spectrum.) 
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FIGURE 13. DOPPLER SPECTRA FOR ACTIVE PASS (PASS 12) X-BAND DATA. 
(Dashed line indicates land spectrum.) 
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orientation, and show Doppler shifts in the same direction of -21.4 
Hz and -18.8 Hz, respectively, for the water spectra relative to the 
land. These shifts indicate horizontal currents of 0.46 m/sec and 
0.38 m/sec in the cross-track direction out of the Straits of 
Georgia. Pass 12 was flown at a different orientation and shows a 
positive Doppler shift of 47.1 Hz, which implies a horizontal current 
of 0.93 m/sec in the same direction. The direction of the current 
indicated by all three passes is consistent with surface observa- 
tions, although the magnitude of the current calculated from passes 
10 and 11 is somewhat too small. 

The system effects which cause the asymmetry in these Doppler 
scans represent a potentially large source of error in the current 
measurement, particularly for the higher-bandwidth data. For Seasat 
data, which has a bandwidth on the order of 1 kHz, the potential 
error could be expected to be quite large. In order to address this 
problem, a set of procedures was developed for making corrections 
for system effects. 

The starting point for this set of procedures was the development 
of a capability for digitally recording the frequency-pl ane measure- 
ments, so that the required manipulations of the data could be more 
easily carried out by computer. This capability was developed as a 
modification of the ERIM image dissector facility, which is normally 
used for digitizing the light distribution in the image plane of the 
optical processor. In order to scan in the azimuthal direction of 
the frequency plane, the Fourier transform lens was mounted on a 
movable platform which was translated by means of a digitally- 
controlled stepping motor. Scanning is also performed electronically 
within the image dissector tube in the range direction, but normally 
the range samples were averaged together. 

Once this digital recording facility was developed, a sequence 
of measurements was made in order to determine the system parameters. 
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which normally vary from one data set to the next. First, the opti- 
cal noise introduced by the processor itself was measured by scanning 
with no film in the signal plane. This set of values was subtracted 
point-by-point (i.e., for each frequency sample) from all of the 
subsequent measurements. Next, a scan was performed with a section 
of signal film which was recorded with no signal input. This data 
was used in two ways: first, it was taken as a measure of the total 
system noise and subtracted from the actual Doppler frequency scans 
for each image subset. In addition, the measurement was used to es- 
timate the modulation transfer function (MTF) of the SAR system. It 
may be assumed that the noise input to the receiver is nearly 
frequency-independent, or at least has a bandwidth larger than the 
SAR system. Therefore, the frequency distribution of the output 
noise spectrum is directly proportional to the frequency response, 
or MTF, of the system. Thus, the noise spectrum was used crs a nor- 
malization factor for the actual Doppler scans in order to compensate 
for the falloff in the system response with increasing frequency. 

During the initial setup before these measurements were made, it 
was discovered that the azimuth telescope introduced an excessive 
amount of distortion into the' frequency plane data. To eliminate 
this distortion, the azimuth telescope was removed from the optical 
bench. This 'necessitated the use of two separate image apertures, 
one in the range focal plane and the other in the azimuth focal 
plane. After the first trial run, pronounced effects due to a non- 
uniform illumination of the signal film were also noted, and steps 
were taken to broaden the illumination pattern. After this adjust- 
ment, the illumination at the edges of the 60 mm aperture was 2.2 dB 
down from the maximum center illumination. 

The data reduction process consisted of four steps. In the first 
step, the data were read from tape and the range samples were aver- 
aged together to yield a single intensity reading for each azimuth 
sample. The resulting set of (typically 3200) samples was written 
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into a disc file for future access. In the second step, the file 
representing total system noise was subtracted from each of the data 
files, and the processor noise was subtracted from the MTF file. In 
the third step, the file representing the system MTF was smoothed 
and truncated at a specified fraction of the maximum value (to avoid 
dividing by zero or a very small number in the next step). In the 
last step, each data file was normalized by the system MTF and 
plotted. Additionally, the last step allows the option of smoothing 
the data and/ or calculating the "median frequency" for each data 
set. The "median frequency," f , is defined as the azimuth fre- 
quency for which the area under the spectrum from f - Af/2 to 
f equals the area from f to f + af/2. The frequency inter- 
val Af was generally chosen as approximately the 3 dB width of the 
spectrum. 

As a test of this procedure, a set of Doppler scans was obtained 
from the 27 July 1978 Porlier Pass data set. Doppler scans were ob- 
tained for an area inside Porlier Pass at an incidence angle of 32° 
and for an adjacent land area. These scans were corrected for the 
system MTF, which was obtained from a section of the signal film 
which was recorded while the transmitter was turned off. The system 
MTF is shown in Figure 14, and the corrected Doppler spectra for land 
and water are shown in Figure 15. The center frequencies calculated 
for these spectra are 237 Hz and 278 Hz, respectively. This 41 Hz 
shift implies a radial velocity of 0.66 m/sec, a horizontal across- 
track current of 1.24 m/sec, and a total current of 1.75 m/sec if 
the direction is assumed to be 45° with respect to the across-track 
direction. This result is in somewhat better agreement with the 
surface observations than the result obtained earlier without the 
full MTF correction. 

The procedure was repeated using the L-band data, resulting in 
the spectra shown in Figure 16. The center frequencies for these 
spectra indicate a small and probably not statistically significant 
Doppler shift of 7.2 Hz. 
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IRCRAFT X-BAND MODULATION TRANSFER FUNCTION FOR CV-580 DATA. 
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FIGURE 15. DIGITALLY RECORDED AND CORRECTED DOPPLER SPECTRA FOR PORLIER PASS 
X-BAND AIRCRAFT DATA. 

(Dashed line indicates land spectrum.) 





FIGURE 16. DIGITALLY RECORDED AND CORRECTED DOPPLER SPECTRA FOR PORLIER PASS 
L-BAND AIRCRAFT DATA. 

(Dashed line indicates land spectrum.) 
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Next, a set of measurements was attempted for the Seasat data 
sets listed in Table 2. The Seasat MTF was measured from a section 
of signal film recorded with a "white noise" input signal, and is 
shown in Figure 17. It may be noted that the system bandwidth as 
indicated by this plot is narrower than the expected Doppler band- 
width of 1.2 kHz. Thus, the full Doppler spectrum was not measurable 
for the Seasat data. It is believed that the system bandwidth is 
limited primarily by the recording system rather than the processor. 

Three Doppler spectra were measured for the Columbia River data 
set (Orbit 150), including one spectrum over the river itself and 
two spectra over the land on either side of the river. The centroids 
of the two land spectra were averaged together and subtracted from 
the centroid of the water spectrum to yield an apparent Doppler shift 
of 9.8 Hz. The water spectrum and the average of the two land spec- 
tra are shown in Figures 18(a) and (b), respectively. The positive 
Doppler shift of the water relative to the land indicates a current 
toward the ground track of the satellite (i.e., west, or out of the 
river) in agreement with surface truth data. The magnitude of the 
shift indicates a radial velocity of 

= I Af = 1 .15 m/sec (23) 

and a horizontal velocity of 

V. = = 3.37 m/sec (24) 

^ sin 20“ 

which is also in fair agreement with the surface truth data. The 
statistical significance of this result is questionable in view of 
the large standard deviation of the land spectra (approximately 26 
Hz). However, the apparent ageement of the calculated and observed 
currents is encouraging. Probable causes of the large variance of 
the land spectra are discussed in Section 3.5. 
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FIGURE 17. SEASAT SYSTEM MODULATION TRANS FE 
THE ERIM OPTICAL PROCESSOR. 
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FIGURE 18(b). DOPPLER SPECTRUM FOR LAND AREA ADJACENT TO COLUMBIA RIVER 
(SEASAT ORBIT 150). 


, The next data set processed was taken from Seasat Orbit 762 over 

the English Channel. Two spectra were measured over water and one I 

i 

over land north of the channel (Figures 19(a)-(c)). The center fre- 
i quencies of the two water spectra relative to the land spectrum were 

measured to be 10 Hz and -14 Hz. These Doppler shifts would imply 
’ horizontal currents of 3.4 m/sec toward the west, and 4.8 m/sec 

toward the east, respectively. The first measurement is in the right 

direction but its magnitude is at least a factor of two too large, 
as compared with the surface truth data. The second measurement is ! 

apparently in the wrong direction. The probable cause of error in 
the second measurement is the highly variable radar return over the 
1 English Channel. As discussed in Section 3.5, such variations in 

. image brightness can cause false spectrum shifts if the signal film ; 

i illumination is not uniform. 

The third set of Seasat measurements was made with data from 
Orbit 651 over the Gulf Stream off the coast of Florida, at a lati- 
tude of approximately 28°N. One spectrum was measured over water 
near the Florida shoreline and a second was taken from an area 
assumed to be inside the Gulf Stream, as judged by the apparent posi- 
tion of the Gulf Stream boundary in the image. These spectra are 
shown in Figures 20(a) and (b). The spectrum from inside the Gulf 
Stream was found to have a center frequency of 4.2 Hz smaller than 
that for the stationary water near the shoreline. This would indi- 
cate a current having a radial component of 0.5 m/sec, or a hori- 
zontal component of 1.4 m/sec in a direction perpendicular to (and 
away from) the ground track. Since the Gulf Stream boundary appears 
to be at about a 45“ angle from the ground track, this would imply a 
total current magnitude of 

l » . i . J . / . s g£ ^ 2.0 m/sec (25) 

sin 45“ 

which is in good agreement with the surface truth information. 
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FIGURE 19(a). DOPPLER SPECTRUM FOR DOVER STRAITS AREA 1 (SEASAT ORBIT 762). 
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FIGURE 19(b). DOPPLER SPECTRUM FOR DOVER STRAITS AREA 2 (SEASAT ORBIT 762). 
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FIGURE 19(c). DOPPLER SPECTRUM FOR LAND NORTH OF DOVER STRAITS (SEASAT ORBIT 762). 
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FIGURE 20(a). DOPPLER SPECTRUM FOR STATIONARY WATER AREA OFF FLORIDA COAST 
(SEASAT ORBIT 651). 



(SEASAT ORBIT 651). 
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All of the preceding Doppler measurements were made on the ERIM 
optical processor as described at the beginning of this section. As 
a test of the feasibility of all-digital processing, a digitally 
processed complex (I and Q) data set was obtained from R. Gray of 
the Canadian Communications Research Centre. This data set consisted 
of an approximately 16 x 16 km segment of Seasat data from Orbit 681, 
processed at full resolution. The image covers the Straits of Juan 
de Fuca and includes a small portion of Cape Flattery, Washington. 
A distinct set of internal waves appears to be propagating into the 
Straits and a gravity wave field is also evident in the image. Test 
areas were selected over land, a uniformly high signal area over 
water, a low signal area over water, and a segment of the internal 
wave pattern. The first three areas consisted of 256 azimuth samples 
and a variable number of range samples. A digital Fourier transform 
was taken for each set of 256 azimuth samples and these spectra were 
averaged together in the range direction to yield the spectra shown 
in Figures 21(a)-(c). The means and standard deviations of the in- 
dividual spectra for the three test areas are shown in Table 3. From 
these statistics, we can estimate the centroids for the test areas 
to be -9.1 ± 0.9 Hz for land, -13.6 =*= 0.6 Hz for the high-signal 

water area, and -10.1 ± 0,8 Hz for the low-signal water area. Note 

2 

that the size ot these test areas are on the order of 1 km as 

2 

compared with the 9 km areas used in the optical Seasat process- 
ing, yet the standard deviation of the centroid estimate is much 
lower. 

The difference between the centroids for the two water test areas 
may be due to a variation in the current or to a systematic error in 
the spectrum measurement which is related to the signal level. The 
latter error could be explained by the presence of an additive "white 
noise" component in the signal, or any type of additive noise with 
zero mean frequency. This would cause the measured centroid to 
approach zero as the signal becomes smaller. Assuming this to be 
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FIGURE 21(a). DOPPLER SPECTRUM FOR LAND AREA IN DIGITALLY PROCESSED 
DATA SET FOR STRAITS OF JUAN DE FUCA (SEASAT ORBIT 681) 
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FIGURE 21(c). DOPPLER SPECTRUM FOR LOW-SIGNAL WATER AREA IN DIGITALLY 
PROCESSED DATA SET FOR STRAITS OF JUAN DE FUCA (SEASAT 
ORBIT 681). 
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TABLE 3 

DOPPLER SPECTRUM STATISTICS FOR THREE TEST AREAS 
IN SEASAT ORBIT 681 DIGITAL DATA SET 


Amplitude 



No. of Range 

Centroid (Hz) 

(arbitrary units) 

Test Area 

Samples 

Mean 

Std. Dev. 

Mean 

Std. Dev. 

Land 

151 

-9.07 

10.98 

1.9 X 10-3 

1.9 X 10-4 

Water 

(higti signal) 

301 

-13.63 

10.52 

2.0 X 10-3 

1.5 X 10-4 

Water 

(low signal) 

151 

-10.08 

10.15 

6.7 X 10-4 

5.7 X 10-5 
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the case, the centroid of the first water test area would be more 
directly comparable to that of the land area, since the mean signal 
levels are nearly the same for these two areas. The Doppler shift 
between these areas is -4.5 Hz, which would imply a horizontal cur- 
rent of 1.5 m/sec into the Straits. This figure is approximately 
three times higher than the tide table prediction for this area. 
Averaging the two water test areas yields a Doppler shift of 2.8 Hz 
and a current of 0.9 m/sec. 

A factor which has not been taken into account in these calcula- 
tions is the phase speed of the Bragg waves responsible for radar 
scattering from the water surface. At L-band, the wavelength of 
these Bragg waves for a 20° incidence angle is 

X = Q -ii 35 . m .. ^ Q^34 m (26) 

^ 2 sin 20“ 

and the phase velocity of these waves is 

Cw =Y^ = 0.73 m/sec. (27) 

Thus, the Doppler shift due to the phase velocity of the Bragg waves 
is 

2C^ sin 20“ 
af^ = — o; 235 ¥“ = 

Assuming the direction of propagation of these waves is in the same 
direction as the large-scale gravity waves visible in the image, 
i.e., into the Straits, this Doppler shift should be subtracted from 
the measured shift before the current is calculated. This results 
in a current estimate of 0.8 m/sec using the high-signal area or 0.2 
m/sec using the average of the two test areas. These values bracket 
the tide table prediction for this area. 

Next, a test area covering the internal wave pattern was exam- 
ined. This area consisted of 1024 x 1024 samples. A Doppler spec- 
trum was calculated for each set of 1024 azimuth samples, and the 
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resulting centroid values were smoothed in the range direction using 
a 100-sarnple boxcar filter. A plot of these smoothed Doppler cen- 
troid values versus the range sample number is shown in Figure 22. 
The signal amplitude is shown for the same set of range sample num- 
bers in Figure 23. Figure 22 shows variations in the Doppler fre- 
quency of about =*=3 Hz correlated with the signal amplitude for the 
first few cycles of the internal wave (on the left of the plot) and 
variations of about ±6 Hz for the two cycles on the right. These 
Doppler shifts would imply a variation in the horizontal velocity of 
±1 m/sec for the waves on the left side and ±2 m/sec for the waves 
on the right side. Comparing the signal amplitude variations in 
Figure 23 with those in Table 3, it is probable that these velocities 
are overestimated by at least a factor of two because of the noise 
effects mentioned earlier. However, a detailed analysis of these 
effects has not yet been carried out. 

No surface truth data is available for these internal waves. 
However, ijr situ measurements on a set of smaller internal waves in 
the Gulf of Georgia indicate velocity fluctuations of ±0.25 m/sec 
i’Gower and Hughes, 1979]. Therefore, the results obtained from the 
Doppler analysis appear to be on the proper order of magnitude if 
the systematic errors are taken into account, 

3.4 ORBITAL VELOCITY MEASUREMENTS 

Theoretical calculations were presented in Section 3.1 which pre- 
dict a broadening of the Doppler spectrum due to gravity wave orbital 
velocities. Figure 24 shows a plot of the amount of this broadening, 
i.e., the spectrum width divided by the width for a stationary sur- 
face, versus the scaled amplitude 


where n is the wave frequency in radians/sec, a is the wave amplitude 
in meters, 0 is the antenna beamwidth and V is the platform velocity. 
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FIGURE 23. SIGNAL AMPLITUDE VERSUS RANGE SAMPLE NUMBER FOR 
INTERNAL WAVE PATTERN IN STRAITS OF JUAN DE FUCA. 
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FIGURE 24. DOPPLER BROADENING FACTOR VERSUS SCALED WAVE AMPLITUDE 
(fv) aircraft X-BAND aircraft data WITH ANTENNA 
BEAMWIDTH OF l.r AND PLATFORM VELOCITY OF 75 m/sec. 
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This plot was generated for the same set of parameters as in Figure 3 
of Section 3.1, although the results are not strongly dependent on 
these parameters. 

Among the Doppler spectra presented in Section 3.4, the only ones 
showing a significant amount of broadening are the X-band Marineland 
spectra. The broadening for these cases range from 1.3 to 2.2 for 
the water spectra relative to the land spectra. Referring to Figure 
24, these would imply a scaled wave amplitude ranging from 0.42 to 
0.84. 

Surface truth data for the Marineland test site indicate a domi- 
nant wavelength of 80 meters and a significant wave height of 1.56 
meters. The frequency of these waves, from linear dispersion theory, 
is 

0.88 rad/sec 

Using a wave amplitude (a) of one half the significant wave height, 
a beamwidth e = 0.019 rad and a platform velocity V = 80 m/sec, the 
scaled amplitude calculated from the surface truth data is 

A = If . 0.45 

This would cause a Doppler broadening of approximately 1.4, which is 
within the range of the observed broadening factors. However, the 
observation of a much larger broadening in some of the cases (Pass 1 
and Pass 3) would seem to imply that the theory presented in Sec- 
tion 3.1 is incomplete, or there is an additional broadening mecha- 
nism that has not been taken into account. Further experimental and 
theoretical work are needed to resolve this question satisfactori ly . 

3.5 ERROR ANALYSIS 

There appear to be at least four factors influencing the apparent 
Doppler spectra of objects imaged by a SAR in addition to the actual 
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motion of the objects. If these factors are neglected or if their 
effects are not distinguishable from actual motion effects, they 
represent potential sources of error in the determination of ocean 
currents and orbital velocities from SAR data. These factors are: 

1. Variations in backscatter cross section with view angle, 

2. Recording system non-uniformities, 

3. Effects of recording and/or processing system MTF and noise, 

4. Aperture and illumination effects in the processing system. 

Each of these factors is discussed in more detail below. 



3.5.1 ANISOTROPIC SCATTERING 

The apparent Doppler spectrum of an object or area will be dis- 
torted if the radar cross section varies over the duration of the 
signal history, i.e., if the cross section is a function of the view- 
ing angle. An evaluation of this effect can be carried out by in- 
cluding the cross section in the expressions developed in Sec- 
tion 3.1. Re-defining A(d) in Eq. (3) as 

A(0) = g(^) (29) 


where a(0) is the radar cross section of the target and g(d) is the 
antenna gain pattern, the Doppler spectrum for a stationary target 
may be written as 



If the angular variation of the cross section is relatively slow, 
a(^) can be expanded in a Taylor series, and 


v'a(O) 1 + ^ ) 


(31) 


For small the above integral can then be evaluated to yield 
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|F(»)| = ^Tcy 


F„(a) - j 


d'(0)v 

2»(0)R„ 


f;(«) 


(32) 


where 1 Fq(o»)| is the Doppler spectrum for a target with unit cross 
section. Assuming a Gaussian form for the antenna gain pattern, as 
in Section 3.1, we can write 


and 



(33) 




o'(0)V 

4d(0)R^b 



(34) 


since b » (FV/sR)^, 


The centroid of this spectrum is located at 


Au) 





8a(0T 


ke^v 


(35) 


Expressed as a fraction of the Doppler bandwidth e' = 2k$V, the shift 
caused by a non-isotropic scatterer is 


Ao) a ' (0 )s 1 Aa 
3' " l6a(0)‘ = 16 7(Ty 


where ao is tne change in the radar cross section during the time 
that the object is in the field of view of the antenna. The degree 
of anisotropy required to produce a shift on the same order of magni- 
tude as that seen in Seasat data is quite certain to occur over land 
targets, and is possible over water as well. The anisotropy of the 
ocean backscatter coefficient depends on the angular dependence of 
the wave height spectrum at short wavelengths, which is not well 
known. Assuming the wave height spectrum to have a cos^ (o/2) de- 
pendence [Tyler, et al., 1974], where 0 is the direction of the Bragg 
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wave relative to the wind direction, the rate of change of the radar 
cross section (including Bragg waves travelling in both the positive 
and negative range directions) is 



-2 sin 9 _cos 0 

sin^ (f) ^ 



(37) 


Thus, for a view angle of 45“ relative to the wind direction, 


g'(e) 

o(») 


0.667 


and the apparent Doppler shift caused by this anisotropy for Seasat 
data (s = 1.73“) is 


^ = 0.0125 

P 


i.e., about 1.5 Hz, which is equivalent to the shift caused by a 
horizontal current of about 0.5 m/sec. 


3.5.2 RECORDING SYSTEM NON-UNIFORMITIES 

An effect similar to that described above could also be caused 
by changes in recording system parameters over the duration of the 
Doppler signal history. Such changes can be produced by variations 
in receiver gain, CRT bias, and/or film response. Film response 
variations can be induced either by emulsion non-uniformities or by 
the development process. Variations are not likely to be more than 
a few percent over a typical signal history, but may be of the same 
order of magnitude as the spectrum changes induced in Seasat data by 
ocean currents. 

The phenomenon of small signal suppression represents another 
possible error source, in which the apparent Doppler spectrum is in- 
fluenced by brightness variations in the scene. The signal film 
density (or transmission) at any point actually represents the super- 
position of many Doppler histories. As long as this superposition 
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is linear, the individual Doppler histories contribute equally {in 
proportion to their scattering cross sections) to the total. How- 
ever, because of film response nonlinearities, the actual effect of 
a bright object's Doppler history is to reduce the effective contri- 
butions of the overlapping histories. This effect will cause the 
frequency spectra for nearby objects to be skewed, resulting in an 
erroneous Doppler shift. Although the effect is more prevalent over 
land, because of the larger range of reflectivities, it may be also 
present to some degree over water. 

3.5.3 RECORDING SYSTEM MTF AND NOISE 

Another factor influencing the apparent Doppler spectrum is the 
recording system frequency response and noise. The film recording 
medium obviously has an upper limit on the spatial frequency that 
can be recorded, and in fact the recording efficiency falls off con- 
tinuously with frequency. This falloff is specified in terms of the 
modulation transfer function (MTF). In order to properly correct 
for this system characteristic, the measured spectrum should be 
divided by the MTF. 

If the system bandwidth is greater than the signal bandwidth 
(i.e., the Doppler spectrum width), a correction can be made for the 
system MTF as outlined above, and the entire spectrum can be used to 
measure the Doppler shift and hence estimate the ocean current. If 
the system bandwidth is smaller than the signal bandwidth, as it may 
be in the case of Seasat data, the MTF correction can still be made, 
but only a portion of the Doppler spectrum can be recovered with 
adecnate signal-to-noise ratio. Thus, the location of the centroid 
and the amount of the Doppler shift is difficult to measure. 

There is also an additive correction factor which should be 
applied to the measured Doppler spectrum to remove the effect of the 
system noise background. Assuming that the noise entering the 
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recording system is frequency-independent, the shape of the noise 
component recorded will be the same as the shape of the MTF. If the 
signal is very large, the noise component wy he negligible, but when 
the signal is low (as it frequently is over water), the noise spec- 
trum may be an appreciable or even a dominant component of the re- 
corded Doppler spectrum. 

3.5.4 APERTURE AND ILLUMINATION EFFECTS 

The azimuth frequency scan is obtained for an appropriate subset 
of the image by setting up the optical processor in a conventional 
manner up to the image plane. Instead of placing a film at tne image 
plane, however, a small aperture is inserted so as to pass only those 
signals which correspond to the subset of the image under considera- 
tion. These signals are passed through an additional set of optics 
which essentially take the Fourier transform of the signals in the 
image plane. The light intensity in the resulting frequency plane 
is then scanned in the azimuth direction and recorded. 

Ideally, the result of this process would be the average power 
spectrum (i.e., the amplitude of the Fourier transform) of the 
Doppler histories of all the objects within the image subset 
selected. The results should not be influenced by objects outside 
•this area, even though their Doppler histories overlap those of the 
targeted objects, and the results should not depend upon the posi- 
tions of the objects within the subset. In order to realize these 
objectives, the following conditions must be met: 

1. The length of signal film illuminated must be at least equal 
to the length of an individual Doppler history plus the width 
of the image plane aperture (in equivalent units), and must 
be centered about the Doppler history of the object at the 
center of the image aperture, 

2. The illumination must be uniform over this area, and 

3. The aperture must be located precisely at the image plane. 
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If the first condition is not met, the high (positive or nega- 
tive) frequency components of some of the Doppler histories will be 
truncated and the frequency distribution will be skewed, resulting 
in an apparent but false Doppler shift. 

If the second condition is not met, the frequency distribution 
of each individual scatterer will be skewed, and the contributions 
from all the scatterers will be unequally weighted. This will again 
result in a false Doppler shift which is dependent on the brightness 
distribution in the scene. 

If the third condition is not met, a perfect discrimination be- 
tween objects inside and outside the image subset will not be 
achieved, and the measured frequency spectrum will be skewed by 
bright objects near the edges of the aperture. It is conceivable 
that this effect may always be present to a small extent {depending 
on the aperture size) due to diffraction effects at the edges of the 
aperture. 

These conditions are all, in principle, controllable and there- 
fore do not constitute a fundamental limitation of the technique 
(with the possible exception of the last-mentioned diffraction 
effect). It may be difficult to control these factors in practice, 
however, and failure to meet the conditions exactly can lead to large 
errors. This is particularly true for Seasat data which typically 
has signal histories occupying a much longer film distance than air- 
craft data. The following set of measurements serves to illustrate 
and quantify this error source. 

The optical processor was arranged for Seasat processing and the 
signal film illumination pattern was widened as much as possible 
using the conventional setup. This resulted in a falloff of 2.2 dB 
in the illumination at the edges of the 60 mm aperture. A signal 
film from Seasat Orbit 982 over Goldstone Antenna was placed in the 
processor so that the image of the antenna was centered in the aper- 
ture, and the Doppler spectrum was measured. Next, the signal film 
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was moved so that the image appeared near the left edge of the 3 kin 
aperture and the spectrum was re-measured. Finally, this was re- 
peated with the Goldstone antenna image near the right edge of the 
aperture. The three recorded spectra (with corrections for system 
MTF and noise) are shown in Figure 25. There is a shift of approxi- 
mately 100 Hz for each of the off-axis spectra as compared with the 
centered spectrum. 

A theoretical explanation of this effect can be developed by ex- 
tending the analysis presented in Section 3.1 to include the effects 
of the signal film illumination function, and to consider signal 
histories located off the axis of the processor. The analysis in 
Section 3.1 assumed that the object was broadside to the antenna at 
t = 0. For an object located an along-track distance x away from 
this point, the signal history is shifted in time by an amount x/V: 
i .e. , 



S^(x, t) 



If the signal-film illumination function is designated as W(t), the 
output of the processor may be written as 


S^lx, t') = J S^(x, t)W(t)q(t - t') dt 

.00 

which has the Fourier transform 

00 

S2(x, u) = Q(o)) J 


( 39 : 


S'! (x, t)W(t) 0 "^“^ dt 


(40) 


where Q(u) is the Fourier transform of the matched filter q(t) 
Assuming the illumination function is Gaussian with width T^, i.e.. 


2 2 
-4t^/Tf 

W(t) = e ® 


(41 


65 


Frequency (Hz) 


FIGURE 25. DOPPLER SPECTRA MEASURED FOR GOLDSTONE ANTENNA WITH 
IMAGE OF ANTENNA (a) AT LEFT EDGE OF APERTURE, (b) 
CENTERED IN APERTURE, AND (c) AT RIGHT EDGE OF APERTURE 
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the above integral can be evaluated to yield the result 



where T = sR/V and = 2kV^. Thus, the amplitude of the 
Doppler spectrum for a point target at x is 

jS 2 (x, u))[ = |Q(u))l exp (w + 



This equation shows that for a point target whose Doppler history is 
centered in the processor (i.e., at x = 0), there is no shift in the 
spectrum although the spectrum is narrowed by the factor (1 + 
T^/T^)~^^^. For a point target located off-axis, there is 
a false shift 



(44) 


in the spectrum. The factor Q(oj) represents the frequency response 
of the processor, which can also cause a reduction in the observed 
Doppler shift. However, this factor can be corrected for by the 
methods discussed in Section 3.5.3 above. The shift caused by the 
aperture or illumination effort is not removable by post-processing, 
although it can be minimized by choosing test areas which have a 
uniform or symmetric brightness distribution. For a test area which 
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is, say, 50 percent brighter on one side of the aperture than the 
other, the spectrum can be approximated by the appropriately weighted 
superposition of the spectra for two point objects at x = =^x /4, 

a 

where x, is the total width of the image plane aperture. The cen- 

a 

troid of this spectrum would be 



tnat is, the error in the centroid divided by the bandwidth 2bT is 


Au) ^a 
2bT " TOVT 



(46) 


If the signal plane aperture width is approximately the same as the 
signal history (i.e., « T) and the image aperture width is 3 
km, as in the case of the Seasat optical measurements reported here, 
this would represent a spectrum shift of about 1 percent of the spec- 
trum width, i.e., 10 Hz. 

Such a shift can easily mask the actual Doppler shift due to 
ocean currents, and probably represents the dominant error source in 
the Seasat measurements. Note that this effect may also be present 
in digitally processed data, depending on the number of azimuth sam- 
ples used in the correlation process. 
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4 

PHASE VELOCITY MEASUREMENTS 

Since a SAR image is generated from a Doppler history which is 
acquired over a period of time on the order of a few seconds, it has 
been proposed that by considering only a portion of each Doppler 
history, a set of images could be generated which correspond to dif- 
ferent time intervals during the acquisition of the data. Moving 
objects should then appear in different relative locations within 
these sub-images. This technique would not apply to fully coherent 
targets, since those having velocities large enough to cause a mea- 
sureable displacement would have Doppler spectra outside of the sys- 
tem bandwidth. However, in the case of gravity waves on the ocean 
surface, the individual scatterers move with the orbital velocity of 
the wave while the positions of the wave crests move with the phase 
velocity, which can be orders of magnitude larger than the orbital 
velocity. 

The theory for this method of measuring the phase velocity of a 
gravity wave is presented in Section 4.1, and a set of measurements 
supporting the theory and demonstrating the technique is presented 
in Section 4.2. 



4.1 THEORY 


The rate of change in the phase of the SAR signal reflected from 
a point object is given by 


•"d 


“ dt " dt 


(47) 


where k is the radar wave number and R is the instantaneous range 
distance to the object. This rate of change is referred to as the 
Doppler frequency of the SAR signal. For an object having a nominal 
range Rq (at t = 0), the instantaneous range is 
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R(t) = Yr^ + s Rg + 




(48) 


where V is the platform velocity. Thus, the instantaneous Doppler 
frequency is 


f 


d 


. “<i . -2V^ ^ 

^tT XRq 


(49) 


If we apply a bandpass filter (i.e., place an aperture in the fre- 
quency plane of the optical processor) to the SAR data such that 
frequencies between f^ and f^ are passed, this is equivalent to 
including only those signals which were collected during the time 
interval and t^ relative to the zero-Doppler time for each 

object in the scene, where 


-xR 


ti = 


2 \l 


^f 

2 


(50) 


If we make an image of the scene using this frequency interval, the 
positions of the objects in the image will correspond to their posi- 
tions in the scene at the time 


- - "i 

‘i = ‘i ^ r <5') 

Where is an appropriate average over the time interval (t^ , 
t^) corresponding to the frequency interval (f ^ , f^) selected, 
and is the along-track coordinate of the object under consider- 
ation. If we make a second image using a different frequency inter- 
val, such that the average of the corresponding time interval is 

Ty, the position of the Uh object in this image will be shifted 
relative to its position in the first image by an amount 

Ari = v.(t2 - t^) (52) 

where v^ is the velocity of the object. 
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For the case of a gravity wave field being imaged by a SAR, the 
moving "objects" are actually collections of scattering centers which 
are themselves moving slowly compared to the phase velocity of the 
gravity wave, but which are varying in brightness in such a way that 
the locus of points having the maximum brightness moves with the 
phase velocity of the wave. This phase velocity is frequently large 
enough to cause a measurable change in the position of the wave crest 
during the imaging time. Thus, the phase velocity ought to be mea- 
surable by determining the shift in the position of the wave crests 
in images generated from different Doppler bands. 



4.2 MEASUREMENTS 


An experiment was carried out to test the technique outlined 
above using L-band SAR data collected by the ERIM X-L system over 
Marineland, Florida, on 14 December 1975. For the particular data 
set (Pass 2) chosen, a gravity wave field with a dominant wavelength 
of 80 meters was present, traveling in the range direction. From 
linear dispersion theory, the phase velocity of these waves is 


“ Vl?- "-2 


(53) 


where \ is the wavelength and g is the acceleration of gravity. 


Four digital images were generated from this data set using four 
different Doppler frequency bands. The frequency intervals used for 
each image are shown in Table 4. A typical Doppler spectrum for this 
data set is shown in Figure 26. This spectrum may be approximately 
represented by the function 


A(f) = e 




( 54 ) 


where f^ = 3.5 hz and 3 = 35 Hz. The average frequency, weighted 
by the signal intensity over the interval (f ^ , f^) is therefore 
given by 
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TABLE 4 

FREQUENCY BANDS AND CORRESPONDING TIMES FOR FOUR 
MAR INELAND IMAGE FILES 


Image 


^2 

T 

t 

File 

(Hz) 

(Hz) 

(Hz) 

(sec) 

1 

0 

11.28 

5.60 

0. 556 

2 

0 

-11.28 

-5.48 

-0.544 

3 

5.64 

16.92 

11.15 

1.107 

4 

-5.64 

-16.92 

-11.03 

-1.095 
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FIGURE 26. DOPPLER SPECTRUM FOR MARINELAND (PASS 2) L-BAND DATA. 
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The average frequencies, and the corresponding average times were 
calculated for each image file from Eq. (50) for a range Rq = 6 km 
and an aircraft speed V = 83 m/sec. The results are shown in 
Table 4. The effective time interval between Files 1 and 2 is 1.10 
sec, and between Files 3 and 4 is 2.20 sec. 


The shift in the position of the waves between Files 1 and 2, 
and between Files 3 and 4, was measured as follows. First, the two 
images in each pair were registered in the azimuthal direction by 
computing the cross-correlation function for the images over a land 
area. This analysis indicated that the maximum correlation was ob- 
tained when File 2 was shifted eight lines (in the azimuthal direc- 
tion) relative to File 1, and when File 4 was shifted 14 lines rela- 
tive to File 3. Next, a set of 200 corresponding azimuthal lines 
was selected from each file. In each subset, five azimuthal lines 
were averaged together (along the wave crests) to improve the signal- 
to-noise ratio, yielding a set of 40 smoothed azimuthal lines. A 
set of 256 range samples was selected from each smoothed azimuthal 
line and the cro^“ correlation function was computed in the range 
direction . 

The cross-correl ation function was calculated both by direct 
computation and by Fourier transform methods. In the direct computa- 
tion method, the correlation function is calculated as 

N / 

p(i) = 2^ S^(j)S2(j i) / 
j=l / 
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where i and j represent range sample numbers, is the signal 
level in the first file, $2 is the signal level in the second 
file, and o-j and 02 represent the standard deviations of the 
signals in Files 1 and 2, respectively. Using Fourier transform 
methods, the correlation function is given by 

p(i) = F.T.-’’ {s*(j)S2(j)} ^ 0^02 (57) 

where S, is the complex conjugate of the Fourier transform of 
Sp $2 is the Fourier transform of S^, and F.T. indicates 
the inverse Fourier transform of the product of these two. In both 
cases, the mean signal is removed from S-j and $2 before the com- 
putation of the correlation function. 

The correlation function calculated from a single line of data 
is an unreliable indication of the image shift because of noise in 
the data. The use of the Fourier transform method for calculating 
the correlation function allows some simple manipulations to be done 
to improve the quality of the results. Since it is the shift in the 
position of the dominant wave which is to be measured, it is desira- 
ble to suppress both the low-frequency signal variations and the 
high-frequency noise, leaving a range of wavelengths around the dom- 
inant wavelength. This is easily done by setting the Fourier coeffi- 
cients to zero outside of this wavelength range. 

The pixel spacing for this data set is 3 meters in range, which 
corresponds to a horizontal distance of 4.04 meters for an incidence 
angle of 48°. Thus, the dominant wavelength is approximately 20 
pixels for the section of data considered (range samples 400-655). 
In the results presented here, the Fourier coefficients corresponding 
to wavelengths greater than 10 pixels and less than 40 pixels were 
set to zero. The correlation function was calculated for each of 
the 40 azimuth lines and averaged over this set. Figure 27 shows 
the resulting correlation function for image Files 1 and 2. There 
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is a clear maximum at a shift of 2.5 pixels, which corresponds to a 
horizontal shift of 10.1 meters. Dividing this by the effective time 
difference of 1.1 seconds; one obtains a phase velocity of 9.2 m/sec. 
The correlation function for Files 3 and 4 (Figure 23) shows a maxi- 
mum at six pixels which corresponds to a horizontal shift of 24.2 
meters. Dividing this by 2.2 seconds yields a phase velocity of 11.0 
m/sec. The latter measurement agrees well with the expected result 
for an 80-meter wavelength, while the former one is slightly smaller 
than the expected result. It should be noted that the zero crossings 
of the correlation function indicate a slightly smaller wavelength 
than 80 meters. The zero crossings occur at a displacement of \/4 = 
18 meters, which would imply a wavelength of about 72 meters and a 
deep-water phase velocity of 10.6 m/sec. In addition, shoaling of 
the waves may cause the actual phase velocity to be smaller than that 
calculated from the deep-water dispersion relation. 
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5 

CONCLUSIONS AND RECOMMENDATIONS 


Evidence has been presented in this report which demonstrates 
that synthetic aperture radar can potentially measure oceanic cur- 
rents, orbital velocities and phase velocities. The evaluation to 
date using conventionally configured SAR systems has shown that there 
are significant errors in these measurements. However, none of these 
systems have been optimally configured for the purpose of measuring 
currents. Desirable system operating parameters for future SAR data 
collection efforts for the measurement of ocean currents should in- 
clude consideration of: 

1 . Incidence Angle 

From purely geometric considerations, it is desirable- to have 
the incidence angle as large as possible in order to maximize 
the radial component of the velocity. However, the radar cross 
section becomes smaller as the incidence angle increases, so the 
geometric considerations may have to be balanced against signal/ 
noise ratio considerations, particularly if the water surface is 
relatively smooth. 

2 • S/W Pl atfor m V elocity 

Since the Doppler spectrum width is directly proportional to 
the SAR platform velocity, it is desirable to keep the velocity 
as low as possible. 

3. A1 titude 

The altitude does not directly affect the width of the 
Doppler spectrum or the amount of the shift, so this is not a 
critical parameter. Stability of the platform is important, and 
the altitude for aircraft SAR's may indirectly affect that param- 
eter, however. 

4. W avelength 

The Doppler shift for a given radial velocity is inversely 
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proportional to the radar wavelength. Therefore, it is desirable 
to operate at shorter wavelengths to maximize the Doppler shift. 

5 . SAR Platfor m Heading 

It would be useful to have at least two passes in orthogonal 
directions in order to resolve the current direction. One of 
the flight paths should be perpendicular to the current to maxi- 
mize the radial component of the velocity. More than two passes 
would be desirable to reduce errors by averaging or to have a 
backup in case of data quality problems. 

6 . Antenna Gain Pattern 

The use of image spatial resolution as the primary design 
criterion for a SAR system leads to a maximization of the antenna 
beanwidth. However, this limits the accuracy of the current 
measurement because the Doppler spectrum width is directly pro- 
portional to the antenna beamwidth. Therefore, if current mea- 
surement accuracy is included as a criterion in the design of a 
SAR system, a tradeoff is required. A1 ternati vely, it may be 
possible to improve the current measurement accuracy while main- 
taining spatial resolution by introducing some structure into 
the antenna gain pattern. For example, a sharp "notch" could be 
created in the antenna pattern by introducing a phase delay over 
one half of the antenna. Also, the tradeoff between spatial 
resolution and current measurement accuracy may be improved by 
the use of advanced systems such as the spotlight technique 
discussed below. 

Both aircraft and satellite SAR current measurements have been 
presented. The majority of the examples presented have utilized 
analog data. Basically, a SAR signal film has been used as the data 
source and lenses are employed to Fourier transform the data in order 
to extract the Doppler information. Recently, SAR digitally cor- 
related data that retains phase (I and Q) information has been used 
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to evaluate SAR current detection. Digital data appears to be advan- 
tageous over analog data in that system effects are more easily 
correctable. 

The measurement of surface currents by the Doppler analysis 
method described in this report is inherently harder using satellite 
data as opposed to aircraft data. This is primarily due to the 
higher velocity of satellite platforms, which causes the Doppler 
spectrum to be wider and hence the measurement of a given Doppler 
shift to be more difficult. In the case of Seasat, this problem is 
compounded by the longer wavelength and smaller incidence angle, 
which cause the Doppler shift for a given current to be smaller by 
an order of magnitude than that for X-band aircraft data. Thus, a 
current which would cause a Doppler shift equal to the spectrum width 
in an aircraft X-band SAR causes a shift of less than one percent of 
the spectrum width for Seasat. It should be noted, however, that a 
satellite platform is more stable than an aircraft platform and a 
one percent Doppler shift may indeed be measurable. 

In view of the advantages offered by satellite platforms, it is 
recommended that studies be carried out to design a more optimal SAR 
system for measuring ocean currents from space. These studies should 
include an optimization of "conventional" SAR systems by the appro- 
priate choice of the antenna pattern, wavelength, and incidence angle 
as discussed previously in this report. They should also include 
consideration of new SAR concepts including bistatic and angle- 
diversity techniques [Walker, 1980] for reducing the Doppler spectrum 
width and resolving the components of the current vector. 

Angle-diversity SAR systems, utilizing the so-called spotlight 
technique, could provide radial 1 ine-of-sight measurements in two or 
more directions in order to resolve the components of the current 
vector. The spotlight technique uses a rotating antenna which "fol- 
lows" a given target area on the surface, rather than a fixed antenna 
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orientation as in a conventional SAR. This allows the system to 
achieve a fine spatial resolution with a relatively narrow antenna 
beamwidth by collecting and processing data coherently over the re- 
quired integration time. However, by choosing to process a shorter 
segment of the data, one can, in effect, narrow the Doppler spectrum 
and therefore increase the accuracy of the Doppler shift measurement. 
Thus, the system provides the flexibility of trading off spatial 
resolution for Doppler measurement accuracy during processing. In 
addition, it may be possible to devise special processing techniques 
which would allow accurate Doppler measurements to be made over an 
extended integration time, and thus preserve a high spatial resolu- 
tion in the current measurement. 

Bistatic SAR systems, utilizing separate transmit and receive 
platforms could also be used to narrow the Doppler spectrum, since 
the spectrum width in this case is proportional to the vector sum of 
the platform velocities. For example, two platforms moving in 
opposite directions produce a Doppler spectrum equivalent to that 
for a monostatic SAR with a velocity equal to the difference between 
the bistatic platform speeds. In both of these cases, tradeoffs must 
be made between spatial resolution and current velocity resolution, 
and a design study would be required to determine the constraints on 
the various system parameters. 

In addition to investigating these new techniques, it is recom- 
mended that a review be made of previously developed SAR moving- 
target indication (MTI) techniques for their potential application 
to the problem of mapping ocean currents. These techniques include 
the use of delay lines for comparing subsequent pulses from the same 
antenna or from different antennas. The ocean current problem is 
somewhat different from the MTI problem in the sense that the motion 
of an extended surface area is involved, rather than the detection 
of a small moving target in a stationary background. Nevertheless, 
some of the techniques developed for MTI systems may be applicable 
to the problem of measuring ocean surface currents. 


82 





RADAR AND OPTICS DIVISION 


2?rii 


This report has also indicated the possibility of measuring grav- 
ity wave orbital and phase velocities under certain conditions. The 
evaluation of these methods has been limited by the small number of 
available aircraft data sets containing suitable gravity wave fields. 
The orbital velocity measurement using the Doppler broadening effect 
has potential applications to the measurement of wave heights, which 
is a topic of considerable oceanographic interest. A further evalua- 
tion of this technique and a study of its possible application to 
satellite data, as is the case of the current measurement, would 
appear to be warranted. The phase velocity measurement is of inter- 
est in verifying wave imaging models but has a smaller range of 
applications and is probably not suitable for use with satellite 
data . 
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ABSTRACT 

Synthetic aperture radar (SAR) has potential 
for mapping ocean surface currents. This instrument 
responds primarily to backscatter from capillary 
waves, which, in conventional SAR processing, are 
assumed stationary. However, these scatterers are 
not stationary; they move with a velocity which is 
the resultant of their own phase velocity plus 
velocities due to the presence of currents and 
longer gravity waves. The radial (line of sight) 
component of this resultant velocity produces a 
Doppler shift in the temporal frequency of the 
return signal, which translates to a spatial 
frequency shift recorded on SAR signal film. The 
relative contributions to this shift by currents, 
gravity wave orbital motions, and capillary phase 
velocities are being studied by a theoretical model 
now under development. 

The authors, using X- and L-band SAR data of 
near shore and Gulf Stream ocean surfaces, measured 
the Doppler shift of moving ocean scatterers rela- 
tive to stationary scatterers. Currents deduced 
from these Doppler shift calculations (averaged 
over a 500 x 500 m^ area) were found to be consistent 
with available sea truth gathered during the Marine- 
land Experiment. Additionally, SAR-SEASAT satellite 
data of the Columbia River, Oregon is being evaluated 
to assess the potential of using SAR to map ocean 
surface currents. 


* The ERIM work reported herein is supported by NOAA Contract No. AOl-78-00-4822 . 
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INTRODUCTION 


This paper discusses the potential use of Synthetic Aperture Radar (SAR) to 
measure surface currents on the ocean. The surface current velocity is measured 
by exploitation of the SAR Doppler signal history, Radial (line of sight) 
velocities of currents shift the Doppler signal history and it is this shift that 
is measured and exploited to obtain the horizontal current velocity. 

Synthetic Aperture Radar is a specialized form of Side-Looking Airborne 
Radar (SLAR) , it utilizes the Doppler history (change of phase) associated with 
the motion of the aircraft, recording both the phase and the amplitude of the 
bacUscattered energy, thus improving the along-track or azimuth resolution. The 
longer the radar data length, the greater the improvement in the along-track 
resolution of the system, 

Synthetic Aperture Radar is a coherent airborne or spaceborne radar that 
uses the motion of a moderately broad physical antenna beam to synthesize a very 
narrow beam thus provldU.g fine azimuthal (along-track) resolution (Brown and 
Porcello, 1969; Harger, 1970), Fine range (cross-track) resolution is achieved 
by transmitting either very short pulses or longer coded pulses which are 
compressed by matched-filtering techniques into equivalent short pulses. 

Usually, the coded pulse is a waveform linearly modulated in frequency, 

Typically, the phase history of a scattering point in the scene is recorded 
on photographic film as an anamorphic (astigmatic) Fresnel zone plate. The para- 
meters of the zone plate are set in the azimuth direction by the Doppler frequen- 
cies produced by the relative motion between the sensor and the point scatterer, 
and in the range direction by the structure of the transmitted pulses. The film 
image is a collection of superimposed zone plates representing the collection of 
point scatterers in the scene, This film is used by a coherent optical processor 
which focuses the anamorphic zone plates into the points which produced the micro- 
wave scatter of the scene (Kozma, et al, 1972). 

The radial motion of an ocean current imaged by a SAR will produce an 
apparent tilt to the phase history as recorded on the signal film. In addition, 
the scatterer history will alsu shift across the signal film. Since the 
synthetic aperture technique utilizes the small range rates of stationary 
reflecting objects to separate their images in the azimuth coordinate, targets 
moving in the range direction appear with altered azimuth positions in the 
processed image, but are still focused, A non-moving target will produce zero 
Doppler when it is on a line perpendicular to the aircraft track. However, if 
the reflecting object is moving toward the aircraft, zero Doppler will be 
produced after the aircraft has passes this perpendicular line, shifting the 
apparent position of the object in the same direction as the flight direction. 
Conversely, if the object is moving away from the aircraft, the apparent position 
will be shifted opposite the flight direction. To illustrate the displacement 
effects which can result from target motion in the range direction, observe 
Figure 1 which shows a SAR image of a moving railroad train displaced from its 
track. This displacement occurs because the radial component of its velocity 
produced radar echo phase rates that are different than the phase rate from the 
fixed targets. The magnitude of the azimuth displacement can be expressed as; 


AX - ^ R_ (1) 

''a ® 


where 


AX 

VR 


displacement , 

line-of-sight velocity component of reflecting object, 
velocity of the radar aircraft, 
radar range of the reflecting object. 
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Note that image displacement is directly proportional to object radial velocity 
and radar range, It is inversely proportional to aircraft velocity and inde- 
pendent of radar wavelength. 

The principle in imaging any surface with a radar is that the backscattered 
microwave energy (echo) received by the radar receiver contains information on 
the roughness characteristics (shapes, dimensions, and orientations) of the 
reflecting area. The principal reflectivity mechanism of a SAR imaging an ocean 
area is via the capillary wave structure. The SAR reflectivity estimate provides 
an apparent brightness that approaches the Bragg-Rice-Philips upper bound 
(Shuchman, et al, 1978), It is therefore necessary to have capillary wave 
structure present before one can attempt to measure current velocities, 
Additionally, the radar observes coherently the projection of 7c “ Vcuj-rents + 
Vcanillaries + ^vave orbital motions onto the slant range vector. The projection 
or Vc is Vr, the coherently sensible range velocity which Doppler modulates the 
SAR signal (Raney and Shuchman, 1978). 

2 . THEORY 

The Doppler frequency shift for a moving target relative to a stationary 
target in a SAR system is 


Af 


D 



( 2 ) 


where V = radial component of target velocity, and 
= transmitted radar wavelength. 

This temporal frequency shift will produce an azimuth spatial frequency shift of 


Af" = 


Afo? 


'AC 


(3) 


on the SAR signal film, 
whefre P = azimuth packing factor, and 
Vac = aircraft velocity, 

Equations 2 and 3 can be combined to relate radial target veclocity to 
Doppler spectrum shift; 


V 


r 


Af"AV 


AC 


TP 


(4) 


This relationship can be used to measure the average radial velocity component 
of an ocean-wave scattering field relative to a fixed-land scattering field. A 
shift in the azimuth spatial frequency spectrum between fixed-land and moving- 
ocean surfaces yields an estimate for Af". Variations in spectrum locations due 
to antenna pointing may be eliminated by choosing imagery having land and ocean 
Imaged simultaneously. 

3 . EXPERIMENTAL TECHNIQUE 

The optical processing configuration used to make the measurements is shown 
in Figure 2. A scanning photomultiplier was used at the output of the optical 
processor to record the azimuth spectrum profile shown in the figure. Note in 
the figure that an image-plane aperture corresponding to an ocean area of 1 km x 



1 km was used in this experiment, Measurements were also made using apertures 
(ocean areas) of 0.5 x 0,5 and 0,25 x 0.25 km. These measurements indicated 
that, although the azimuth Doppler scans become noisy as smaller ocean areas are 
used, acceptable scans can be obtained using an aperture corresponding to 0,25 x 
0.25 km in radar space. Additionally, the SEASAT measurements to be discussed 
were made with a 3 x 3 km aperture. 

Figure 3 presents the X-band (3.2 cm) azimuth profiles for 1 km x 1 km 
land and ocean areas imaged on Pass 3 at Marineland on 14 December 1975 (Shemdin, 
et al, 1978). Both areas were centered 6.8 km from the end of the pass, The 
land area was near nadir and the sea area was just inside the 60° incidence angle 
bounce line. As the figure indicates, the Doppler shift is to the left with 
respect to the land, thus indicating that the radial motion sensed is toward the 
radar (see bottom of Figure 2) . The mid-point of a line across the 3 dB points 
of each Spectrum was plotted for each scan. The separation of the mid-points was 
approximately 1.55 divisions. This separation indicated a shift of -8 line-pairs/ 
mm (1 div. = 5.2 dp/mm) between the two scans. Using this value for Af" in Eq , 

(4) yields 


_ 8 X 103(0.032) (76 m/sec) 
2P 2 (13,327) 


0,75 m/sec 


(5) 


Figure 4 presents the X-band azimuth profile for Pass 7 (180° heading with 
respect to Pass 3). Note that the spectrum has shifted an amount approximately 
equal but in an opposite direction to that of Pass 3, 

Figure 5 is an L-band (23.5 cm) example of a near-shore scan showing that 
L-band can also be used to detect radial motion. The shift distance is not as 
great as for X-band. However, similar velocities can be measured, Recall that 
Eq . (4) is a function of radar wavelength (X) and azimuth scaling factor (P) . 

It should be noted that, since the shift for a given velocity is larger on the 
azimuthal Doppler scans for X-band than for L-band, the X-band scans are 
potentially more accurate, due to the larger scale in comparing stationary and 
moving Doppler scans. 


4 . EXPERIMENTAL RESULTS 

A series of motion scans of near-shore surface currents were made at L-band 
and X-band. Additionally, motion scans of the Gulf Stream were also made at 
X-band only. 

Figure 6 is a graph of the L-band near-shore measurements per pass vs. the 
spectrum shift from the DC reference (see Figure 5). In the figure, the letter 
"L" represents land or stationary scans while S and D represent wave scans made 
from data collected in shallow and deep water, respectively. Although this data 
indicated spectrum shifts from stationary to moving targets, the data was disap- 
pointing, in that neither the land measurements (L) nor the ocean data clustered. 
This rather disappointing result was felt to be caused by the operation of ERIM's 
L-band radar at Marineland* rather than an inherent limitation of SAR in sensing 
currents. 

Alternatively, Figure 7 presents X-band spectrum shifts as a function of 
the three passes. The X-band results were encouraging, in the sense that 
stationary and moving Doppler scans clustered and Passes 3 and 7 (180° with 
respect to each other) shifted in opposite directions, Note that Pass 1 showed 
no azimuthal shift induced by velocity. This confirms the limited surface- 
current truth provided by CERC (Army Corps of Engineers) . The CERC data 
indicated that an 85 cm/ sec current existed on 14 December at Marineland, with 
a direction of travel toward shore. 


* For a description of the ERIM ,X-L SAR system see Rawson, et al (1975). 



X-band data from the Gulf Stream was also evaluated, Scans A-H (locations 
shown in Figure 8) were made using a 1 x 1 km aperture, Figure 9 compares 
azimuth scans of location A and location G. A graph plotting the location of the 
scans on the Gulf Stream data versus Af from the DC reference is shown in Figure 
10. If the average Af between the Gulf Stream and non-Gulf Stream data is 
calculated using Eq, (4), a radial velocity (V^) of 0.63 m/sec results. The 
direction of this velocity is determined by examination of Figures 8 and 9, 
Unfortunately, gravity-wave components were present during data collection. 

These waves were traveling in the same general direction as the Gulf Stream. 

Thus, there is some question as to whether the velocity sensed by the SAR could 
be an orbital velocity wave component. The authors do not believe this to be the 
case, however, because deep-water waves were measured both in the Gulf Stream and 
out of it; thus, any orbital motion effect from the gravity waves was the same on 
each side of the Gulf Stream boundary. 

5. SEASAT CONSIDERATIONS 

NASA launched in June 1978, an oceanographic satellite called SEASAT. One 
of the five sensors on board the satellite was a SAR. This SAR, which imaged the 
ocean at a range of 850 km produced imagery 100 km wide by up to 4000 km in 
length. The incident angle of the SEASAT SAR is nominally 20°. The resolution 
obtainable is 0 m azimuth by 25 m ground range. 

SEASAT data collected between June and October 1978* is being evaluated by 
the authors in respect to its ability to measure surface currents. Data from 
Revolution 150 (July 7, 1978) is being evaluated in detail. In particular. 
Revolution 150 imaged the Columbia River, Oregon. Figure 11 shows the resulting 
SEASAT imagery. The image represents subswath 2, its resolution is on the order 
of 10 m azimuth, 25 m ground range. The scale factor is nominally 1:215,000, thus 
the image represents an area of —52 x 35 km. 

This image offers the advantage of attempting to calculate the velocity of 
the Columbia River (reported to be 2. 5-3.0 m/sec at the time of flight) by either 
the use of Eq. (1) or the Doppler centroid method. A 3 m/sec velocity using 
Eq. (1) would result in a displacement of 115 meters. This corresponds to shift 
of approximately 0.5 mm on tne image, a distance too small to accurately discern. 

The Doppler method was also used in an attempt to measure the velocity of 
the Columbia River. Four land scans were compared to four Columbia River scans 
(different areas) . The four data values resulting from the four pairs of 
measurements produced a horizontal velocity of 4.8 m/sec. This value is somewhat 
higher than the reported 2. 5-3.0 m/sec velocity. All SEASAT measurements of the 
currents area to date appear to have higher than reported velocities. The cause 
for this apparent bias is not well understood at present. 

6 . SUMMARY 

The above described technique utilizes the measurement of Doppler shifts in 
the SAR signal history induced by radial velocities. Thus oceanic surface 
currents to be measured must be traveling perpendicular to the SAR flight track. 
Additionally, because a radial velocity is sensed (line-of-sight) the measure- 
ments must be corrected to the horizontal datum plane. This is done knowing the 
SAR geometry. 

Azimuth traveling currents can theoretically be measured by observing 
defocusing effects in the SAR signal processor (Shuchman and Zelenka, 1978). 
However, the authors believe the defocusing effects are too subtle to be 
accurately measured. 

Surface currents are inherently harder to measure using satellite data as 
opposed to aircraft. This is due to the increased Doppler bandwidth of the 
satellite in respect to aircraft. For example, the aircraft Doppler bandwidth 
at Marineland is approximately 100 Hz, while SEASAT 's is approximately 1000 Hz. 


* On October 10, 1978 SEASAT suffered a catastrophic power loss. 



A 1.5 m/sec radial velocity produces a Doppler shift of approximately 10 Hr thus 
a shift measurement for the aircraft is a 107» change while for the satellite the 
change is 17«. It should be noted, however, chat the satellite platform is more 
stable than an aircraft's and a 1” measurement may indeed be measurable. 

The effect of a gravity wave field on a current area to be measured has not 
been studied in detail to date. Initial observations indicate the wave motion 
causes a syminetrical broadening of the Doppler history and therefore does not 
alter the shift induced by the current. 

In conclusion, the current measurements presented herein appear consistent, 
in the sense that Doppler shifts were detected and the shifts appear in the 
proper direction, However, the limited sea truth available for Marineland 
prevents a definitive statement as Co the ultimate feasibility of using SAR to 
sense current motion. Additionally, high velocity current data has been 
collected off Vancouver Island, which is a physical situation very different 
from the large area, low velocity situation of Marineland, This data will be 
analyzed in the near future. The SEASAT measurements need further theoretical 
consideration to explain the overly large Doppler shift in the Columbia River 
bar example. 
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1 , ABSTRACT 

Synthetic aperture radar Doppler histories have been used to 
obtain estimates of ocean surface currents. For currents in excess 
of 1 meter/second, airborne radar measurements agree with surface 
observations to within 25 percent. Satellite measurements, however, 
yield current speeds far in excess of surface estimates. The tech- 
nique developed to obtain the Doppler information is described, and 
possible sources of error are discussed. 

2 . INTRODUCTION 

This paper discusses the potential use of synthetic aperture 
radar (SAR) to measure ocean surface currents. SAR systems are 
characterized by an azimuth (along-track) resolution comparable to 
that in the range (cross-track) coordinate. The azimuth resolution 
is achieved by using the Doppler history associated With the plat- 
form motion to synthesize an antenna (aperture) many times longer 
than the physical antenna; this in turn produces a correspondingly 
reduced effective radar beam width. In effect, the Doppler infor- 
mation is employed to locate targets in the azimuth coordinate. 
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The Doppler shift is positive as the platform approaches the target, 
zero as the target is directly abeam, and negative as the target 
recedes. 


However, such processing assumes stationary targets, and in the 
case of ocean SAR imagery this assumption does not hold true. Rather, 
the primary radar scatterers are short gravity and capillary waves 
whose net velocity is the vector sum of their own phase velocity, 
the spatially and temporally periodic orbital velocities of longer 
gravity waves on which they may be borne, and the larger scale, more 
uniform velocity of any surface currents which may be present 
(Shuchman et al, 1978). The radial component of this resultant ve- 
locity vector (i.e. its projection along the radar line of sight) 
will Induce an additional Doppler shift in the radar return signal. 

In principal, this additional shift should be detectable by compar- 
ison with the Doppler history of stationary targets. 


3 . THEORY 

SAR Doppler frequency shifts for moving targets relative to 
stationary targets are given by Shuchman et al (1979) 


Af 


D 



( 1 ) 


where 


V = radial component of target velocity 
A*' * transmitted radar wavelength 

This temporal frequency shift will produce in turn an azimuth spatial 
frequency shift on the SAR signal film given hv 


Af ' 


where 




( 2 ) 


V 

P = — = the azimuth packing factor 

V = platform speed 

v^ = recording film speed 

Combining Equations (1) and (2) yields an expression which can be 
used to e$tlmate the average radial velocity component of an oceanic 
scattering field 


- Af *AV 
'^r " 2P 


(3) 
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The required estimate of Af’ can be obtained by the experimental 
technique described in the next section. 

In conventional SAR prooesslpg, the target induced Doppler shift 
causes a corresponding shift of the apparent azlmutih position of ^he 
target in the image; thus » target movement toward the radar will 
Induce positive Doppler, and the target will appear to be farther 
ahead of the platform than it really is. Direct measurement of 
these azimuth displacements may yield estimates of surface currents, 
especially in regions of very high shear (Shemdin et al, 1980). For 
a target at range R, the magnitude of this azimuth displacement is 
given by 


V 

Ax - ^ R (4) 

Note that this displacement is directly proportional to target vel- 
ocity and range, but inversely proportional to platform velocity 
and independent of radar wavelength. 

4. EXPERIMENTAL TECHNIQUE 

The optical processing configuration used to make the Doppler 
measurements is shown in Figure 1. A scanning photomultiplier was 
used at the output of the optical processor to record the azimuth 
spectrum profiles such as those shown In the Figure, The estimate 
for Af in Equation (3) is obtained from the difference in peak 
frequency for scans corresponding to land and ocean imagery. Errors 
in the location of each peak due to antenna pointing can be eliminat- 
ed by comparing land and ocean data obtained simultaneously, i.e. 
having identical azimuth coordinate. Note in the Figure that an 
image-plane aperture corresponding to an ocean area of 1 km x 1 km 
is Indicated. Aircraft measurements were also made using apertures 
(ocean areas) of 1.0 x 1.0, 1.0 » 0.5, 0.5 x 0.5, 0.5 x 0.25 and 
0.25 X 0.25 km2. These measurements Indicated that, although the 
azimuth Doppler scans become noisier as smaller ocean areas are used, 
acceptable aircraft scans can be performed using an aperture corres- 
ponding to 0.25 X 0.25 kni2 in radar space. The SEASAT measurements 
to be discussed were made with a 3 x 3 km^ aperture. 

5. RESULTS 
5.1 Aircraft 

I In July of 1978, several flights were made in the Straits of 
Georgia in the West coast of Canada, by the Canadian CV-580 research 
aircraft equipped with the ERIM dual polarized X- and L-band SAR 
(Rawson et al, 1975). The site selection for the experiment was an 
attempt to maximize the ocean surface current signal and minimize 
complications due to the presence of gravity waves. Narrow channels 
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Fig. 1: Optical processing configuration for azimuth 

Doppler frequency scans. 


in this region, such as Active Pass and Porlier Pass, are protected 
from Pacific Ocean waves by Vancouver Island, and are characterized 
by peak tidal currents commonly in excess of 2 meters /second. 
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AIRCRAFT SAR FLIGHT GEOMETRY 
FOR SURFACE CURRENT ESTIMATES 


Fig. 2: Aircraft flight geometry for SAR -derived current 

estimates. 

Figure 2 shows the SAR aircraft flight geometry for the X-band 
surface current estimates. Figure 3 Is a sketch which Indicates the 
resulting velocity measurements obtained for Porller pass. Note 
that eight scans were made of (stationary) land to compare to (non- 
statlonary) water. All surface current estimates by SAR are compared 
with surface observations in Figure 4. In all cases, the aircraft 
SAR correctly sensed the direction of the current but underestimated 
the magnitude. This bias may be the result of comparing spatially 
averaged SAR estimates with relatively higher point measurements 
obtained In situ near the center of each channel. 

5.2 SEASAT 

SEASAT, an experimental oceanographic satellite, was launched 
In June 1978. One of the five sensors on board was an L-band SAR, 
which Imaged the ocean at a range of 850 km and incident angle of 
about 20°, producing imagery 100 km wide by up to 4000 km In length, 
and an obtainable resolution of 8 m azimuth by 25 m ground range 
(1 look). A SAR image of the Columbia River obtained 7 July 1978, 
and one of the English Channel obtained 19 August 1978 will be ex- 
amined here. 

The Columbia River Image is shown in Figure 5. The resolution 
Is approximately 10 m in azimuth and 25 m in ground range. The scale 
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Fig. 3; SAR-derivecl surface current estimates for Forller 
Pass. (Radial velocities deduced from alrbortle X-band SAR.) 
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Fig. 4: SAR-derlved surface current estimates for Active 

and Porlier Pass compared with in situ observations. 

factor of the original imagery is nominally 1:215,000, and represents 
an area of about 52 x 35 km^. Tide table predictions indicate that 
the current at the mouth of the Columbia River was oriented approx- 
imately perpendicular to the satellite track, with a speed of about 
2.5 - 3.0 m/sec at the time of image acquisition. 

The radial component of a 3 m/sec velocity vector perpendicular 
to the SEASAT track is about 1 m/sec. For this value of v , Equation 
(4) predicts an azimuth displacement of about 115 meters for a moving 
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Imi'th 



Fi^. 5: SF-ASAT SAR inuige of the Columbia River entrance. 


target; this corresponds to a distance of about .5 iran on the original 
Imagery. For a point-like moving target, such as a surface vessel, 
the velocity distribution of the scene I 3 discontinuous and the as- 
sociated displacement is detectable. In contrast, ocean currents 
are typified by smooth, continuous /'locity distributions, and the 
"targets" are patches of ocean water with spatially varying average 
radial velocities. It is unclear at the present time what the final 
radar signature would be as a result of such velocity distributions, 
and more work is needed to develop these ideas as a potential tech- 
nique for surface current estimation. 

The Doppler method was used in an attempt to measure the vel- 
ocity of three sections of the Columbia River (Figure 6). The three 
measured areas produced radial velocities of 1.5, 4, and 5.5 m/sec. 

The latter pair of values, when translated to the horizontal, are 
much higher than the reported 2.5 - 3.0 m/sec velocity. The cause 
of these large overestimates is being investigated. 

The image of the English Channel shown in Figure 7 was acquired 
at a time corresponding to a predicted tidal current with a compon- 
ent perpendicular to the SEASAT track of 0.5 to 1.5 knots. The 
prominent linear features which appear in the Channel apparently 
represent surface manifestations of Channel bottom topography, as 
they correlate well with submerged sand waves which appear on local 
bathymetric charts (Alpers, in prepn.). To investigate the pos- 
sibility that such features in the SAR imagery are induced by the 
surface current velocity distribution, scans were made across the 
Imagery to measure both relative image intensities and Doppler shifts. 
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Fig. 6: SAR-derived surface current estimates for the 

Columbia River. 



Fig. 7: SF.ASAT SAR image of the Thames Estuary and English 

Channel between Clacton and Calais. (Orbit 762) 

The results are shown in Figure 8. There does appear to be signifi- 
cant correlation among the linear bathymetric features (as Indicated 
by the intensity scan) and the surface current distribution (as 
implied by the Doppler shift scan) . 

6. DISCUSSION 

The ocean current measurement techniques we have described here 
yield only a single component of the current. Thus, measurements 
must be made from at least two different viewing angles (preferably 
differing by 90 ) in order to resolve the full current vector. Such 
measurement pairs can be acquired almost synoptically by aircraft, 
whereas a single satellite system would collect such data only over 
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a number of hours or days, and thus might be restricted to the 
measurement of loi|iger tnrni mean currents. In principal, it is pos- 
sible to estimate I a second current component perpendicular to the 
radar line of sight by measuring the degree of defocuslng induced; 
however, these effects may be too subtle for accurate measurement 
(Shuchman and Zelenka, 1978). The application of the hydrodynamic 
equation of continuity to a field of radial current components may 
also lead to resolution of the complete vector field; however, the 
initial radial vector field must probably be of much higher accuracy 
than that presently obtainable by this method. 


The short discussion which follows outlines the more signifi- 
cant sources of experimental error and their relative importance 
with regard to aircraft and satellite applications. 


6.1 System Sensitivity (Geometry and Doppler Bandwidth) 


The radial component of an ocean surface current is given by 
(refer to Figure 2) 


v “V sine.slnG, 
r A i 

where 

V = the magnitude of the current 
6^ = the current direction referred to azimuth 
0^ = the radar incidence angle 

The aircraft data examined here was obtained at incidence angles of 
30 to 60 degrees, while SEASAT viewed the ocean at a less favorable 
angle of about 20 degrees. 

Furthermore, the ocean current signal v and the associated 
Doppler shift Af^ for a moving target referred to land must be 
detected against the background of Doppler shifting Induced by the 
platform motion (V) Itself. Assuming that motion is uniform, the 
bandwidth for a SAR system with azimuth resolution 6 is given approx- 
imately by (Tomiyasu, 1978) 



of the system sensitivity to ocean surface currents 
the deviation Af^ to this Doppler background, given 


• sln0^) V sin6^ = s • v sin6^ (5) 

factor s is shown in Table 1 for the SAR's used. 


Thus, a measure 
is the ratio of 
by 


^ = /-ii 

B ' AV 


The sensitivity 
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Table 1 


The Factor s (Sensitivity to Ocean Surface Currents) 
jFor CV-580 and SEASAT SAR 


SAR 

System 

‘ & 
(m) 

X 

(m) 

! V 
(m/s) 

°i 

(deg) 

s 

(m/s)“^ 

CV-580 X-band 

3 

.032 

125 

45 

1.1 

CV-580 L-band 

^ 3 

.23 

125 

45 

0.15 

SEASAT L-band 

10 

.23 

7000 

20 

0.0042 


By this measure, we see that the CV-580 X-band channel is in- 
herently about 7 times more sensitive than the L-band channel, and 
250 times more sensitive than the SEASAT system. Values of s 
larger than 1 cause image suppression due to bandwidth limitations, 
since this implies that the moving targets are Doppler shifted en- 
tirely out of the scene. This is a problem which can be avoided by 
special processing. 

Note also that, in principle, an X-band satellite system with 
azimuth resolution of aboyt 350 meters would have an s factor equi- 
valent to that of the CV-580 X-band channel. 

6.2 Platform Stability 

The true velocity vector of an aircraft or satellite-borne 
radar antenna is not constant in time. In particular, changes in 
the magnitude and direction of the antenna velocity can result in 
non-zero velocities and accelerations in the radial direction; this 
in turn can cause shifts in the Doppler signal history which con- 
taminate the current-induced shifts being sought. Satellite platforms 
are more stable than aircraft, although it should be noted that some 
aircraft are equipped with an inertial system which may partially 
compensate for these contaminations, 

6.3 Surface Gravity Waves 

Interaction of the short gravity and capillary wave scatterers 
with the almost periodic orbital motions of longer gravity waves 
may broaden the Doppler spectrum and/or make it asymmetric. Pre- 
liminary observations have detected such asymmetries and broadening, 
though a detailed study of these effects has not been completed. 

It should be noted that generally low sea state conditions existed 
during the collection of data summarized in Figure 4. 

6.4 Anisotropic Backscattering 

In the course of synthesizing a SAR aperture, a particular 
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target region is viewed by tlie radar from differing angles. If the 
backscatterlng cross section of that region is anisotropic, the 
associated Doppler spectrum will not be uniformly weighted for all 
frequencies. Pertinent examples of such anisotropic scatterers 
might be surface gravity waves, or the banks of tidal channels. 
Experiments are planned to determine the effects of this anisotropy 
on the radar phase history. 

7 . SUMMARY 


Airborne and SE^SAl’ SAR Doppler histories have been examined 
as a means of estimating ocean surface currents. 

For currents in ; excess of 1 meter/second, airborne X-band SAR 
measurements of the current speed were found to be about 25 percent 
lower than the corresponding surface observations. Thik discrepancy 
may be due in part to the fact that SAR estimates are spatially 
averaged Instantaneous measurements, while the in situ observations 
represent temporal averages taken at a point. Additional errors 
may be Introduced into the SAR estimates by platform instability, 
the presence of surface gravity waves, and target regions with aniso- 
tropic backseat taring properties. 


The SEASAT L-band SAR is poorly suited for this type of 
measurement, and estimates obtained with this system were in sub- 
stantial disagreement with tidal current predictions. The poor 
SEASAT results ma y be due in part to unfavorable viewing geometry, 
large Doppler bandwidth Induced by high platform speed, and the 
effects of gravity waves and anisotropic backscatterers. A future 
satellite SAR system designed for measuring ocean surface currents 
might be improved using a shorter radar wavelength and a coarser 
azimuth resolution. 


While the aircraft 'results reported here are encouraging, much 
more experimental and theoretical work is needed to determine the 
effects on SAR Doppler spectra of such things as long wave orbital 
motions, anisotropic backscatterers, and spatially varying surface 
velocity distributions. 
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1 . ABSTRACT 

The propagation of gravity waves and the wave-current inter- 
action phenomenon has been studied with SEASAT synthetic aperture 
radar (SAR) data obtained along the east coast of the United States 
near Cape HatteraSi North Carolina. Particular attention is focused 
on SEASAT Revolution 974, where the meteorological conditions, wave 
generation In deep water and the resulting wave field are discussed 
and compared to a numerical shallow water wave model. 

The resulting coastal wave data results show a reasonable 
representation of the refraction pattern near Cape Hatteras. Optical 
Fourier Transforms (OFT) and Fast Fourier Transforms (FFT) generated 
In deep and shallow water are compared with the resulting measured 
and modeled wave field. Results indicate that the SAR is a valuable 
instrument for current and wave detection and particularly for wave- 
current interaction studies. Wave length versus water depth inversion 
techniques have shown encouraging results for using SAR data in deter- 
mining bathymetric structures. 

2 . INTRODUCTION 

Remote sensing techniques have recently become an increasingly 
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lmp«*r t .ir*' tiuJ In tlu- otf.inli i . 1 1 1 y in tin* (ii'tfction 

of lu «‘.in w.tvfK, iiirronth .iiul o«f.init bottom ioat iiros . Tlio Inui^fry 
I'ol lot tod from .i spot i.i 1 i /od form of a Hide looking; im.i>;inK radar 
caliod SAR fli'wn i>n SHASAI' was oxaminod for tboHo pfiont>r»*n.i . A 
limitod d.ita sot of SAR lm.i>’t-ry obt.iinod for Rovolution 974 along the 
ea«t roaKt of the I’nlted States, near Cape lintteras. North Carolina, 
showed part iitil.tr ly interesting wave patterns int hiding interaction 
or the refraetion «»f snrta«e waves by the (iulf Stream .ind shoaling 
due to the bathvmetry (Figure 1). 

An.ilysis of OH prodiu ed estinjiites of wave length and wave 
direction for several sites in deep water, in the Culf Stream and in 
shoaling waters. Ihese d.it.i and the FFT information of all dominant 
wave periods and direct ions were input ft>r a wave refraction and 
ocean current numeric.il model. The resulting wave rays, including 
the interaction of the w.ivi- field with the ocean current and the 



Fig. 1: SF.ASAT SAR Digital Iirvage (100 x 100 km area) 

of Cape Hatteras Region. 
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shallow water refraction were compared to SAR FFT's and surface 
measurement. 

3. OPTICAL AND DIGITAL FOURIER TRANSFORMS 

Anticipating the necessity of optically processing SEASAT SAR 
datA, ERIM investigated means of optically extracting wave informa- 
tion. The methods that appear to hold the most promise are based on 
the two-dimensional Fourier transform. Wave parameters which can be 
Inferred from the Fourier transform of the corrected SAR data are 
wavelength, wave direction, and spatial uniformity (straightness of 
the wavefronts, amount of refraction, number of wave directions, and 
continuity of wavefronts). The optical two-dimensional transform is 
simply achieved by Fraunhofer diffraction, as reported by Shuchman, 
et al (1977). 

In addition to approximately 80 OPT's made of the Cape Hatteras 
data, eight fast Fourier transforms were made on the dltlgal data 
(Figure 1) processed by the Jet Propulsion Laboratory. The range 
coordinates of the digital data were analytically corrected for slant- 
to-ground geometry (Feldkamp, 1978). Subsections of 9.2 km x 9.2 km 
with 18 meter resolution were extracted from the digital data. The 
18 meter pixel spacing were converted to 36 meter data by 4 pixel 
averaging in order to Increase the coherence in the image. The 
average value of each azimuthal line was subtracted from the line to 
remove the trend of intensity falloff with Increasing range distance. 
Two-dimensional fast Fourier transforms were performed on each 256 x 
256 cell subsection to yield raw directional wave number spectra with 
a Nyqulst wave number of 0.09 . The raw spectra were smoothed by 

replacing each value with the average of the surrounding 5x5 cell. 
The approximate number of degrees of freedom for the resulting 
spectrhm is 142 (Kinsman, 1965). The 99% confidence limits are then 
il.5 dB (Jenkins and Watts, 1968). Figure 2 is an example of deep 
water spectra obtained approximately 80 km due west of Cape Hatteras. 

4. WAVE REFRACTION 

The Cape Hatteras data (Revolution 974) demonstrated the capa- 
bility of the SEASAT SAR to image shallow water gravity waves under- 
going refraction, and the Interaction of the Gulf Stream with the 
wave field. 

A graphical wave refraction diagram was drawn using the local 
bathymetry of the Cape Hatteras region and compared to the SAR 
derived wave orthogonals (Figure 3). The solid lines on the figure 
were constructed using classical methods outlined by Johnson, et al 
(1948) with the deep water wave length and direction Information 
being obtained from the ERIM OFT's. The dotted lines were SAR 
derived, coming from the OFT and FFT analysis. Care was taken to 
begin the refraction landward of the Gulf Stream boundary so that 
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Fig. 2: Contour Plot of Fast Fourier Transform of 

Deep water Spectra. 


the effect Gulf Stream current on wave refraction would not enter 
Into the comparison. In general, the comparison between the refrac- 
tion diagram and the SAR derived orthogonals Is good. 


To carry the refraction study one step further, a wave refrac- 
tion and ocean current refraction model developed by Hayes (1977, 
1980) was employed to compare the results of the OFT and FFT data to 
surface measurements. 


The basis for the ocean current Interaction study lies In early 
work by Johnson (1948) which pointed out that major ocean currents 
such as the Gulf Stream may have an appreciable effect on the height, 
length, and direction of waves approaching the shore and, under some 
circumstances, may cause almost complete reflections. 


Also, Kenyon (1971) points out that tills pheonomenon of wave- 
current Interaction may cause trapping of waves at certain Incident 
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Fig* 3: Actual SAR-Derlved (dotted lines) versus 

predicted (solid lines) orthogonals for SEASAT 
Revolution 974 (depth in feet on figure). 


angles. He also shows that the radius of curvature is a character- 
istic length scale for the refraction and that it decreases with 
increasing current vorticity and with decreasing group speed relative 
to th^ current. He continues, "For waves that propagate with 
(against) a variable current, the rays will bend in the direction of 
decrejasing (increasing) current speed. Therefore, under the right 
condijclons, it is possible for Waves that propagate with (against) a 
current to be trapped about a local minimum (maximum) in current 
speed. It is also possible that waves that propagate with (against) 
a current could be reflected from a local maximum (minimum) in current 
speed." 

Longuet-Hlgglns and Stewart (1960) show the change in angle of 
a refracted wave to be the same as that reported by Johnson (1948); 
that is , 

sin 9 = sin 9n 

(1 - (V/Co) sin 00 )2 

where 9 is the refracted angle, Bq Is the incident angle, V is the 
current speed, and Cq Is the incident phase velocity. 
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LONGITUDE 

Fig. 4: Sample wave refraction diagram without the 

Gulf Stream. 

OFT and FFT periods and directions were first input to the wave- 
current model and wave rays were produced without consideration for 
refraction due to the Gulf Stream (Flgu_e 4). The resulting angles 
and wave spectra wore linearly recombined (Kinsman, 1965) in shallow 
water and along the ray paths and compared to ottier FFT spectra. 

The results indicated that the magnitude and direction of the model 
runs did not compare favorably to the FFT's produced in the Gulf 
Stream and in shoaling waters. Therefore, the rays were inltlallzecl 
again, this time including a non-uniform current distribution (Hayes 
1980). The resulting rays (Figure 5 and 6) illustrate the bending 
of the wave rays due to the shear produced by the wave-current inter- 
action. When ithe recombined shallow Water spectra, which included 
the modeled Gulf Stream, were compared to the SAR FFT’s in shallow 
water, the directional resolution compared mere favorably than cases 
without the modeled current. 

Therefore, the results illustrated here indicate that the SAR 
apparently has the ability to image the effect of ocean currents as 
well as of shallow water and deep water refraction. 
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Fig. 5: As for Figure 4 with the Gulf Stream included, at 

a period of 10.40 seconds. 
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Fig. 6: As for Figure 4 with the Gulf Stream included, 

at the peak period of the NOAA Data Buoy (11.1 seconds). 
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